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1 Overview 



Since the introduction of electromagnetic wavelets in |K94j . I have been intrigued by two 
related questions: (a) what are the sources responsible for their emission and absorption, and 
(b) can these sources be realized and used to launch and detect the wavelets? In principle, 
it should be easy to find the source for a known solution of a linear system. However, these 
wavelets are pulsed-beam propagators that are, in a 'poetic' sense at least, waves generated 
by a single event in complex spacetime as seen by an observer in real spacetime. Their 
sources, therefore, should be (ignoring polarization for simplicity) extended delta functions 
5(x' — z), where x' is the real observation point and z the complex source point. Attempts 
to make mathematical sense of all this have presented a number of challenges that have 
occupied me for the past several years. I am now pleased to report that the effort has been 
successful and the results are interesing. The sources for scalar ('acoustic') wavelets have 
been rigorously computed both in spacetime and Fourier space and will be presented here. 
The spacetime expressions turn out to be singular objects with a wealth of geometric struc- 
ture and, most interestingly, their Fourier transforms are extremely simple — contrary to 
all my previous expectations. To understand how these simple expressions generate such 
rich spacetime structures, I complete the circle by computing the pulsed beams from their 
Fourier sources. As an unexpected byproduct, this gives an angular spectrum representation 
of time-harmonic complex-source beams, used extensively in the engineering literature, that 
generalizes Hermann Weyl's well known representation of the fundamental solution for the 
Helmholtz equation. The sources for electromagnetic wavelets are polarization and magneti- 
zation densities obtained by multiplying 5(x' — z) by electric and magnetic dipole moments 
or, more generally, convolving them in time with variable dipoles representing a ' driving 
signal.' 

1 believe that question (a) above has now been largely addressed, though some important 
points of interpretation remain, as do no doubt many others. This opens the possibility 
suggested by question (b). Furthermore, although the wavelets were introduced in classical 
electrodynamics, I have hoped from the beginning that they may be useful in QED. This 
is encouraged by the simplicity of their momentum space propagators, which could offer a 
useful computational tool. 

2 Waves, wavelets, and complex spacetime 

I begin with a brief review of my past efforts to extend classical and quantum theories to 
complex spacetime and interpret the results physically. By that I mean that the imaginary 
spacetime coordinates, and any other extras associated with analyticity, are to be understood 
directly in terms of common observable attributes and not merely as a technical device for 
proving theorems or exotic higher dimensions inaccessible to mortals stuck in the 'real' world 
like the poor souls in Plato's cave. 

I have tried not to impose an a priori grand vision but, rather, interpret the imaginary 
coordinates in each theory by understanding their effects within that theory. Consequently, 
the interpretations vary somewhat from one theory to another. But they all have in common 
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the following theme. In the extended theory, certain singular points (evaluation maps on 
Gelds or wave functions, source points, etc.) become 'inflated' to extended objects. This 
transformation is determined by analyticity and the particular theory. In every case, the 
structure of the objects is shaped by the equations of the theory and their degrees of freedom 
are specified precisely by the complex spacetime coordinates. The real coordinates give the 
center, and the imaginary coordinates the extent and orientation of the object in space and 
time. 

These ideas are similar in spirit to wavelet analysis, where a function of one variable (time, 
say) is expressed in terms of an additional variable describing the scale or resolution in the 
first. This analogy goes farther in the treatment of massless than massive fields, since the 
latter have an intrinsic scale and thus cannot be scaled arbitrarily. For relativistic fields 
with mass, spacetime 'orientation' includes velocity, and this makes the complex spacetime 
an extended phase space. The relativistic coherent-state representations for massive Klein- 
Gordon and Dirac fields constructed in |K77|. IK78j (for single particles) and |K87|. [K90J 
(for free quantized fields) interpolate between 'time-frequency' and 'wavelet' descriptions, 
behaving like the former in the nonrelativistic regime and like the latter in the ultrarelativistic 
one. In fact, there is a very close correspondence between the nonrelativistic limit in physics 
and the narrow-band approximation in signal theory; see \K90\ IK94} K96j. 

Although the results cited in this section are not new, I believe they have acquired some cur- 
rency because of substantial progress recently in the understanding of the sources associated 
with retarded holomorphic 1 fields. The new results focus on massless fields, but it is likely 
that similar computations exist for massive fields where the integrals are more difficult. 

Sources describe the breakdown of analyticity due to natural singularities and physically 
necessary branch cuts. What I find especially fascinating is that such branch cuts behave 
much like 'real' matter. Depending on the theory, they carry charge, mass and spin, and they 
emit and absorb radiation. In spite of their simple origins, they turn out to have surprising 
and complex (pardon the expression) properties, the pursuit of which has the feeling of 
exploring hitherto unknown forms of matter and not merely the mathematical properties of 
branch cuts. The results of this search have intrigued and inpired me, and I hope to share 
this excitement with the reader. 

Partial reports have appeared in [KOOl IKoTl iKTTTal IK021 lK02allK03j . but the detailed com- 
putation of complex spacetime point sources and their Fourier transforms, and the angular 
spectrum representation of complex- source beams, has not appeared previously in print or 
preprint form. 

2.1 Spacetime and Fourier notation 

Real spacetime vectors will be written in the form of complex Euclidean vectors with real 
space coordinates and an imaginary time coordinate: 



1 According to most referrences, holomorphic simply means analytic. Possibly the term is used to empha- 
size complex as opposed to real analyticity — or perhaps just to intimidate the reader! 




2 2.2 
x = r —t, 



r = \x\. 




3 



This is known in physics as Hot' (we take c = 1) and is often regarded to be unphysical 
because it cannot be used in 'generic' curved spacetimes where the metric tensor cannot 
be continued analytically in t. (See [MTW73, p. 52], Farewell to "ict.") Therefore I feel 
compelled to explain at the outset why I use it nevertheless. The view of physical time as 
an essentially imaginary variable comes from the idea of complex distance [KOO], rooted in 
Euclidean spacetime, which is the basis for the recent progress in the analysis of sources. In 
previous work I had followed a more conventional path, starting with known quantities in 
Minkowski space and extending them analytically when possible. That approach, reviewed 
later in this section, works best for free fields, where the splitting into positive and negative 
frequency components provides a natural setting for analytic continuation. To understand 
sources, we must rather look at propagators, i.e., waves emitted by a single point source 
S(x), since waves emitted by all other sources can then be obtained by convolution. However, 
propagators are much more complicated than their elliptic counterparts, the potentials due 
to point sources in Euclidean space. In the hyperbolic case, we must choose between retarded 
and advanced propagators (or hybrid ones, like Feynman's). Each choice is singular not only 
at the source but also on at least one half of the light cone, and it is not clear whether and 
how such solutions can be extended to complex spacetime. The Euclidean case, by contrast, 
is a picture of simplicity. Because the 5-source in M n is spherically symmetric, its potential 
G n depends only on the distance r. To extend G n anlytically we need only extend r, and 
this is easily done. The only complication is that the resulting complex distance r(z) is 
double-valued in C n because of the square root in its definition, and a branch cut must be 
chosen. For even n > 4, the extended G n is even in f and so does not depend on the choice 
of branch, while in all other cases it does. 

This suggests an alternative strategy for extending physical fields without excluding the pos- 
sibility of sources. Since the original signature is irrelevant once spacetime is made complex, 
why not begin in the Euclidean setting, where operators are elliptic and life is simple, and 
continue analytically to the Minkowskian world? There may be advantages to starting di- 
rectly with Euclidean concepts, unencumbered by baggage imported from Minkowski space, 
even if in the end we intend to study propagators. Thus we begin with the extended poten- 
tial G^(z) in C 4 and look for extended versions of the retarded and advanced propagators 
D + (x) and D~~(x). However, since the concepts of propagation and causality are foreign to 
the Euclidean world, G±(z) cannot distinguish between retarded and advanced propagators. 
In fact, it turns out to be the extension of the Riemann function, which is the sourceless sum 
of the advanced and retarded propagators (point source). To split off the retarded part, a 
branch must be chosen for the spatial complex distance f 3 , representing the 'distance' from 
the complex source point at z E C 3 to the observer at x' e M 3 . This gives a splitting into 
extended propagators, D + (z) and D~(z) that depend on the branch cut of f 3 . This branch 
cut, T> C M 3 , is a 'blown up' version of the original point source which acts as a disk source 
for D + (z) and D~(z). 2 Finally, D ± (x) are recovered by an 'edge of the wedge' type limit 

D ± (x) = Urn \ D ± (x - iey) - D ± (x + iey)) (2) 

e\0 I J 

2 Going through the cut means entering a world where distance becomes negative and D ± (z) trade places. 
As z becomes real, the disk contracts to a point and the door is closed! 
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with y in the future cone of M 3,1 . This shows that while analytic extension is very simple 
for Euclidean potentials, it is rather involved for propagators, lending support to the idea of 
beginning in the Euclidean setting. But in this approach, physical time must be imaginary, 
just as 'Euclidean time' must be imaginary when beginning with Minkowskian or Lorentzian 
spacetime. 

At any rate, the dismissal of Hct' in |MTW73] may have been premature. Even in general 
relativity, analytic continuations in time and space have borne some rich fruit, even if the 
physical basis of the procedure is often ill-understood. For example, an exquisitely simple 
geometric derivation of the Hawking temperature for Schwartzschild black holes is obtained 
HI79 by analytically continuing the metric in time, interpreting the 'Euclidean time' coor- 
dinate as an angle, and choosing its period to make the horizon a coordinate singularity like 
the origin in polar coordinates. The reciprocal of the imaginary time period is interpreted in 
the usual (KMS) way as a temperature, and this turns out to be nothing but the Hawking 
temperature! I confess that I do not understand this derivation in more than a formal way, 
but analytic continuation has, in any case, become the main strategy of black-hole thermo- 
dynamics, as explained in|Kr03j. The analytic continuation of spatial coordinates also has 
an honorable history in relativity, having played a major role (and conceptually an equally 
obscure one) in the discovery of charged spinning black holes by Newman et al. |N65j ; see 
also |N 731 INW74| IK01a| IN02j . And then there are the theories of twistors and H-spaces (see 
|FS71lHHP'J.7Sl IBFPHOj and references therein). 

Returning to our summary of notations, we denote by M' the 'momentum space' dual to M 
and write the pairing between the two as the complex Euclidean scalar product, 

k-x = k-x — ut, k = (k,iu) e M 1 , x=(x,it)eM. (3) 

The Fourier transform and its inverse are written 

F(k) = [ dx e~ ik - x F(x) (4) 
F(x) = I dk e ik - x F{k), (5) 

JM 1 



where we have introduced the notation 

dx = dx dt dk 



dk du 



d=d/2% (6) 



which eliminates most factors of 2n in Fourier analysis. Strictly speaking, we should include 
factors of i in front the integrals 3 since the time and the frequency are imaginary, but this 
seems overly pedantic and we choose to leave them out, accepting the consequent possibility 
of an extra factor of % now and then. 

Note: Even though x = (x,it) and k = (k,iu), we will write F(x,t) and F(k,u) instead 
of F(x, it) and F(k, iuj) to keep the notation simple. 

3 To be precise, one should begin with dual vectors (fe, — iu>'), since the frequency in units of it is —iu>', 
so that dx = idxdt and ctk = —id~kd~io'. The change to ui = —uj 1 (or fc 4 = — ki), which makes the phase 
velocity +k/u>, does not affect d~k since d~u>' ■ ■ ■ — d~u> ■ ■ ■ . 
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We also need the following cones in M and M' . 
• The positive and negative frequency light cones 



C± = {(k, iu) G M' : ±u = k > } where « = |fc| . (7) 



Note that we exclude the 'DC component' k = 0, i.e., nonvanishng constant solutions; a 
wave must oscillate to be a wave. (This is analogous to the admissibility condition in wavelet 
theory.) 

• The future and past cones 



\y\ 



V± = {(Vi i u ) £ M : ±u > a} where 
which are characterized by the duality relations 

y EV ± ^ k-y <0 Vk eC ± . (9) 

• The solid positive and negative frequency cones 

Vl = {k= (k, iu) E W : ±uj > k > 0} (10) 
which are the convex hulls of C±, characterized by 

k E Vl & k ■ y < G V± . (11) 

• The double cones 

c = c + uc_, v = v + uv_, v' = vluv\ (12) 

Our main complex spacetime domains will the forward and backward tubes 

T± = {z = x - iy : x E M, y E V±}, (13) 

and their union, which we call the causal tube, 

T = {z = x - iy : x E M, y 2 < 0}. 

2.2 The extension of relativistic fields 

All free relativistic fields extend analytically to T in a sense to be explained. This can be 
done directly in spacetime by the analytic-signal transform (AST) |K90[ IKS92| |K94j 4 



p(x-iy) = — I -F{x-sy/u), y=(y,iu)EV. (14) 

2m > /_ 00 s — iu 



4 The AST introduced in the above references differs from (|14|l by a sign in order to be consistent with 
hyperfunction theory. It will be seen to be Lorentz-invariant when applied to free fields. 
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To see how this works, suppose to begin with that F is a general function not necessarily 
satisfying any differential equation, and substitute the Fourier expression for F: 

Fix -iy) = — / / dk < ik -< " Flk). 

2« y_ 00 s-iu J M , 

Assuming the order of integration can be reversed, compute the integral over s by closing 
the contour in the upper or lower half-plane depending on the behavior of the exponential. 
This gives 

00 fit 

e~ isk - y/u = uG(-k ■ y) e k ' y (15) 

where is the Heaviside function and we have introduced a simple notation for the sign 
function, which will be used frequently: 



u = Sgnw and 0(0 = I ' ^ (16) 



(This is easily remembered since u just a ID version of a unit vector u.) That gives the AST 
as an extension of the inverse Fourier transform: 

F(x-iy)=u [ dk Q(-k-y)e iHx - iy) F(k), (17) 

JM' 

which also tells us that the Fourier transform in x of the AST is 

F(k,y)= I dx e~ ik - x F(x-iy) = uQ(-k-y)e k - y F(k). (18) 

J M 

In general F is not analytic, in spite of the notation, because the Heaviside factor spoils 
analyticity as y varies. (F may not have an analytic extension.) However, if F(k) is supported 
in V (fTUj) . then the Heaviside function in (JTTj) disappears after correlating positive frequencies 
with the future cone and negative frequencies with the past cone: 



1, k 6 VI and y G V± 
0, keVl and y G V T 



Q(-k-y) = Q(uju) = 

and we get a third expression for the AST, 

F(z) = ± I dk e ik ' z F(k), zeT ± . (19) 



Now 

k G Vl , y eV ± =^ k ■ y < n(a =p u) < 0, 

hence the factor e k ' y in the Fourier-Laplace kernel e tk ' z decays exponentially. If F(k) does not 
grow exponentially, the integrals (|19j) define holomorphic functions in T ± . (Note, however, 
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that the exponential decay gets weaker and weaker as y approaches the light cone. This will 
later give the ability to focus pulsed-beam wavelets.) 

Of course, it cannot be claimed that F is an extension of F itself since only the positive- 
frequency part is represented in T + and only the negative-frequency part in 71. Define the 
partial boundary values 

F(x =F iO) = lim F(x =F iey) , y E V + 

e\0 



which do not depend on the particular choice of y. Then F is a boundary value of F in the 
same sense as (j2J), 5 

F(x) = F(x - iO) - F{x + zO). (20) 



Generally, the restrictions of F(z) to the disjoint domains T ± are unrelated holomorphic 
functions, but if F(x) vanishes in an open region of spacetime, then ()20j) implies that they 
are part of a single holomorphic function. This is the famous edge of the wedge theorem 

It is instructive to compute the AST of a function depending only on time: 

1 f°° ds n , 
F(t-iu) = — I - — —F(t-s), u^O 

F(s) 



2m 

1 f" f/.s 



2iri ./-oo t — iu — s 



u 



/oo 
duj Q(uu)e- iuj{t - tu) F(uj) (21) 
-oo 



or 



. f dtu e- 1 ^-™) F(u), u>0 
F(t-iu) = t {°n . u . . } h 22 

This is a pair of analytic signals extending the positive and negative frequency parts to the 
lower and upper complex half-planes, a useful concept introduced by Dennis Gabor [G46J 
which also explains the name of our transform. (Actually, Gabor worked with real signals, 
where it suffices to consider only one of the above pair since the other is merely its complex 
conjugate.) 

Equation (|2~Tj) shows that F is a convolution of F with the Cauchy kernel, 



F(t)= ds C(t - s)F(s), T = t-tu } C(t) = , (23) 

J-oo 2mr 

and that the Fourier transform of C is 

POO 

C(lo, u)= dt e iuJt C(t -iu)=u e(cjM)e-" u . (24) 



3 The sum gives a spacetime version of the Hilbert transform [K 
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The role of the sign u in the AST can be understood by assuming that F(t) is compactly 
supported in an interval /. Then, according to the edge of the wedge theorem, F(r) is 
holomorphic for all r ^ / and we may therefore rewrite (J21|) as a contour integral 



1 f ds 



T 



2ixi L, t 



F(s), T = t-iu, (25) 
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where 7 is a closed contour surrounding /, running in the positive direction at s — iO and 
in the negative direction at s + iO. Equation (|2~K|) is Cauchy's formula for the values of an 
analytic function outside of 7 if the contour can be deformed through infinity in the Riemann 
sphere. This interpretation would not exist without the sign in the definition of the AST. 

Equation (|23j) shows that the sign u in (JTTj) gives a positive orientation to the boundary of 
the lower-half complex time plane (i.e., the time axis as seen from below) and a negative 
orientation to the boundary of the upper-half complex time plane (the time axis as seen from 
above), which is obviously correct. This interpretation carries over directly to the spacetime 
AST. In a reference frame where y = (0,u), (fTIj) reduces essentially to (|2"T|) : 

F[X,T) = - — r / F[X,S), T = t — lU. 



2vri 

In the sequel, it will be helpful to keep in mind the following correspondence between the 
4D geometry of spacetime and the ID geometry of time: 

Vl , C± < > {to : ±uj > 0} 

V± < > {u : ±u > 0} 

T ± < — > {t^iu : ±u > 0}. 

2.3 Massive fields and relativistic coherent states 

Free particles and fields of mass m > satisfy the Klein-Gordon equation 

DF(x) = (A — d 2 t )F(x) = m 2 F(x). (26) 

(In the case of spinor fields, all components satisfy (|26j). We consider scalars for simplicity. 
Dirac particles and fields are treated in |K87| lK90j .) In Fourier space, this means 

(k 2 + m 2 )F(k) = 



F(k) = 2vr5(m 2 + k 2 ) f(k) = 2tt5(E 2 - uj 2 ) f(k), E = Vm 2 + k 2 
= I {5(E - u)f(k, E) + 8{E + u)f(k, —E)} 

for some function / defined on the double mass shell 

fi ± = {k G M' : ±uo = E}, £1 = fi + U . 
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Therefore 



F{x) = I dk e ik - x f{k) + I dk e ik - x f(k) = I dk e tk - x f(k), 
Jn+ Jq_ Jo. 



where 



*=2E = «iSE (27) 
is the Lorentz- invariant relativistic measure on Q ± |IZ80j . Applying the AST gives 

F{z)=u [ dk®{-k-y)e lk - z f{k) 
Jn 

= ± I dke ik - z f(k), z = x-iy eT±. 

This has been used to build a coherent-state representation for massive particles (where f(k) 
is a 'c-number' function vanishing on Q_ to give a positive-energy solution) and fields (where 
f(k, E) and f(k, —E) are creation and annihilation operators). These representations have a 
physical interpretation relating y to the expected energy-momentum of the coherent states. I 
am not going to repeat this construction here as it is readily available |K77| IK781 IK871 IK90j , 
but merely indicate how the above interpretation comes about. 

Fix y G V + and consider the exponential e k ' y as a function of k G Q + . The Lorentzian scalar 
product satisfies the reverse Schwartz inequality 



k ■ y < —mX, where A = \/—y 2 = Vu 2 — a 2 > 
which becomes an equality if and only if k is parallel to y: 

k-y = -mX <^> k = (m/X)y = k y (28) 



e 



k - y < e - mX and e k ' y = e~ mX <S> k = k 



y ■ 



Therefore e k ' y acts as a ray filter in momentum space, favoring those plane waves propagat- 
ing approximately in the direction of y. The larger we take A, the stronger the filter and the 
more collimated the ray bundle passed by it. The coherent states are defined in momentum 
space by 



e,(A0 = e-** , (29) 

so that they act as evaluation maps on the Hilbert space 7i of holomorphic solutions with 
inner product defined in L 2 (dk): 

F(z) = {e x \f), (Fx\F 2 ) = (fi\f 2 ). (30) 

This makes Ti. a reproducing kernel Hilbert space, and that kernel is 

K(z',z*) = (e z ,\e z ) = [ dk e ik < z '- z ^ = -iA+(m; z' - z*), (31) 
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where A + is the Wightman 2-point function for the scalar field of mass m. 

To obtain a resolution of unity, choose any oriented 3-dimensional submanifold S in M as a 
prospective configuration space and the 6-dimensional submanifold 

ax = {x - iy e T + : x E S, y 2 = -A 2 }, A > 

as the associated classical phase space, with A > arbitrary. A symplectic form and covariant 
measure on <j\ are chosen as follows. Define the invariant 2-form and 6-form 

a = dx^ A dy^, a 3 = a A a A a. 

Then the following are proved: 

• The restriction of a to <j\ is a symplectic form if and only if S is nowhere timelike, i.e., it 
must be locally spacelike or lightlike. In other words, a\ is a reasonable phase space if and 
only if S is a reasonable configuration space, so the symplectic geometry is compatible wth 
the spacetime geometry. 

• If S satisfies the above condition, then a 3 defines a positive measure du.\ on ax and we 
have a resolution of unity 

dfxx(z) \e z )(e z \ = I, d/j x = C(X,m) a 3 (32) 

where C(X, m) is an invariant and I is the identity operator on L 2 (dk). 

• The physical interpretation of a\ as a phase space is confirmed explicitly in the case when 
S is flat by the expected positions and momenta in the state e z , 

(Pf,) ez =ay^ y li = -1mz fl fi = 0,1,2,3 (33) 

( x j(t))e z = Xj x= ix.it) j = 1,2,3 

where a(A, m) is an invariant and Xj are the Newton- Wigner operators in the Heisenberg 
picture at time t. (When S is curved, the positions operators obtained by quantization on 
ax do not commute; see |K76j .) 

• The parameter 1/ A measures the uncertainty or resolution in the momentum of coherent 
states parameterized by z G a x , in accordance with the above discussion of ray filters. 

• Unlike the usual spacetime representation, the coherent-state representation admits a 
conserved, covariant probability current density, given by 

Uz) = J^ = iHz y a m- iHz) a J^L. (34) 

JM ' dyv K J dx» K J dx» K J 

• In the nonrelativistic limit one obtains a coherent-state representation of the centrally 
extended Galilean group with Gaussian measure in momentum: 

= C e- my /u dxdy, y = —(P). (35) 

m 
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The weight function is imposed on the relativistic (hence uniform) measure dfj,\ in compen- 
sation for the deformation of the mass shell Q + to a 3-plane at infinity (the nonrelativistic 
momentum space). Upon applying a 'holomorphic gauge transformation' jKM80j to the non- 
relativistic wave functions (solutions of Schrodinger's equation), the representation becomes 
identical to the Bargmann- Segal representation of the Weyl-Heisenberg group: 

^nr e mz2/4u F NK rf/i A NR -> Ce- mlzl2/2u dxdy. (36) 
2.4 Massless fields and wavelets 

In the limit m — > 0, the 'reconstruction constant' C(X,m) in (}3*2*j) diverges and the coherent- 
state representation is no longer square-integrable due to disappearance of the 'mass gap' 
k > m > 0. At the same time, the symmetry group grows from the Poincare group to the 
conformal group C. In its realization as SU(2,2), C acts on T ± by matrix- valued Mobius 
transformations. Thus it is reasonable to look for resolutions of unity adapted to the new 
symmetries. This took me several years to realize, and only when studying wavelet theory 
in preparation for Ingrid Daubchies' 1990 Ten Lectures conference |D92] did I understand 
that scaling needed to be brought into the picture. Since massless fields are important 
in classical as well as quantum physics, I decided to begin with classical fields, the prime 
examples of which are acoustic and electromagnetic fields. This led to the construction of 
acoustic and electromagnetic wavelets |K92|. IK941 IK94aj . From a foundational as well as 
applied point of view, I believe the electromagnetic wavelets hold far more promise and 
the 'acoustic' ones serve mainly to simplify the analysis by stripping away all complications 
related to polarization. I now briefly review the construction of acoustic wavelets, leaving the 
electromagnetic ones to the end of the paper where they and their sources will be constructed 
from the scalar acoustic sources. 

Solutions of the scalar wave equation are given by 

□F(x) = F(x) = I dk e ik ' x f{k), (37) 

Jc 

where dk = d~k/2n is the massless version of (|27j) on the light cone C (|12j). Applying the 
AST gives the extension to T±, 

F(z) =u [ d~k B(-k ■ y)e ik - z f(k) = ± [ dk e ihz f{k). (38) 
Jc Jc± 

In the massive case, we chose the coherent states to be the complex conjugates of e lk ' z . 
However, this cannot be done now because e tk ' z does not vanish near uj = and this will 
spoil the resolution of unity. (In wavelet terms, e tk ' z is not admissible.) We get around this 
difficulty by changing the inner product of solutions to 

(F 1 \F 2 ) = (f l \f 2 }= / - fl ( k yf 2 (k), u>0. (39) 

Jc K 
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For v — this is the Lorentz- invariant inner product, but it will turn out that we need v > 1 
to obtain admissible wavelets. Next, write the extension (J38j) in the form of an inner product 

C fju 

F(z)=u - v K v Q(-k-y)e ik -*f(k) = (rl> x \f), (40) 
Jc K 

where 

ipz{k) = uk u <d(—k ■ y) e~ ik ' z * (41) 
are to be the 'acoustic wavelets' in Fourier space, representing spacetime solutions 

\& z ( x >) = u [ dk k u e{-k ■ y) e ik < x '- x+ ^ = ty(x' - z*). (42) 
Jc 

From the invariance of dk under k — > — k it follows that 

= -*(*), (43) 

hence it suffices to compute \1/ in T + . For any fixed y G V + , the function \l/(x — iy) may be 
called a 'mother wavelet' from which all others are obtained by complex translations. Notice 
that imaginary translations generalize scaling. 

As in the massive case, where a great deal of freedom existed to choose a phase space a due 
to the abundance of coherent states, there is now a lot of freedom in choosing a family of 
\l/ 2 s to build a resolution of unity. Perhaps the simplest choice for a continuous frame is by 
analogy with ID wavelets, which are parameterized by position and scale. 'Position' is now 
x, and we take 'scale' to be u, since it dominates the other scale parameters y. Thus we fix 
any time t, say t = 0, and set y = (this will give spherical wavelets). Our parameter space 
is then 

E = {z = (x, -u):xe M 3 , u ^ 0}, (44) 

which is Euclidean spacetime consisting of real space and imaginary time coordinates, with 
the Euclidean time u acting as a scale in M 3 ' . A quick dimensional analysis shows that to 
compensate for the weight k v in Fourier space, we need the measure 

dfx u (z) = C v dx \u\ u ~ 2 du, 
with C v adjusted to give a resolution of the identity l v in L 2 (dk/ ' k v ): 

[ d(i v (z)\* z )(* B \=I v , C v = 2 v /T(y-1). (45) 

J E 

Thus we must take v > 1 to get 'admissible' wavelet representations, and (j43|) then represents 
solutions of the wave equation as superpositions of spherical wavelets centered at x with a 
pulse duration T = 0(u). These wavelets are sourceless, converging onto x when t < and 
diverging from x when t > 0, and their radius at the waist t — is, like the pulse duration, 
R = 0(u). As promised, u controls all scales. 
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Converging spherical wavelets are unnatural under ordinary conditions. 6 We want to elim- 
inate the converging (advanced) part and retain only the diverging (retarded) part. The 
resulting wavelets will have sources, but the splitting cannot be done by brute force (e.g., 
multiplying by the Heaviside function) since that will spoil the analyticity and amount to 
introducing sources with infinite support. A natural separation into advanced and retarded 
wavelets was found in |K94| Chapter 11] while computing *&(z). We review this because it 
foreshadows the recent developments. 

Assuming v is a nonnegative integer, we have for z G T + 



UK 1% ' 

C+ 

ik-. 



#(*)=/ dk k v e lk ' z = (-d u ) v G(z) (46) 

1 



G(z) = dke 



c+ Atc 2 z 2 ' 

where G(z), the original inadmissible kernel with v — 0, is most easily computed by using 
the Lorentz-invariance of the integral. Now 

— z 2 — (t — iu) 2 — (x — iy) 2 = (t — iu — r)(t — iu + f) 

where 

f(x — iy) = yj (x — iy) 2 , Re f > (47) 

is the complex distance from the source point iy to the observation point x, which will be 
studied in detail below. Thus G(z) can be expanded in partial fractions, giving 

g(*) = J-( — 2 — _ 1 — -1 ze r + 

8iT z r yt — iu — r t — tu + r) 
V(z) = ty + (z) +®-(z) 

*±(z) = ± ir{U + 1) - . (48) 

The expression for ^f(z) remains valid if v is not an integer, provided an appropriate branch 
cut is chosen. It splits naturally into retarded and advanced parts without spoiling holomor- 
phy everywhere, as would a brute-force splitting. Instead, the two parts acquire singularities 
from the complex distance function. For given y ^ 0, 

• they diverge on the branch circle S of radius \y\ in the plane y^, where f = 0; 

• they are discontinuous across the branch disk T> spanning S; 

• T> is a source for ^ f± , canceling in 

• and are pulsed beams converging to, then diverging from T>. 

We have called r the 'distance' from an imaginary source point iy to a real observation point 
x. At this stage, such language must be viewed as 'poetry' since the idea of a point source 



6 There are exceptions, including the recent method of time-reversed acoustics where sound is recorded 
by an array of microphones which are then played in reverse, sending the waves back. 



14 



at iy has not been denned and it is not even clear what it means. In physics, a complex 
distance function identical to r was the basis for the construction of spinning, charged black 
holes; see |N65| INT3] . In engineering, 'complex- source pulsed beams' similar to the above 7 
have been applied extensively since the 1980s, and their time-harmonic components, known 
as 'complex-source beams,' since the 1970s; see |HF01j for a comprehensive review. But until 
recently, no serious study seems to have been undertaken to make mathematical sense of the 
idea of a "complex point source," whether time-harmonic or pulsed, and therefore of how such 
beams may be realized. (Sec [HLKOOJ for an early attempt.) Perhaps this is because their 
singular and convoluted structure in the near zone appears to make the requisite analysis 
difficult if not impossible. We will see that the sources are tractable in spacetime and, most 
significantly, even simple and computationally effective in Fourier space. 

Although the requirement v > 1 precludes a Lorentz-invariant wavelet representation for 
acoustic waves, it does admit one for electromagnetic waves, where the invariant measure on 
C has v = 2. Note that 

~, = ^ < 49 > 

is scaling invariant as well as Lorentz-invariant. In fact, the associated Hilbert space of solu- 
tions of Maxwell's equations carries a unitary representation of the full conformal group C, 
as proved by Gross |Gr64j . The electromagnetic wavelet representation is likewise covariant 
under C [K94, Chapter 9], and this opens up some interesting applications. I will not discuss 
the details here since more recent developments are discussed later. 



3 Point sources in complex space 

We now begin implementing the 'Euclidean strategy' of setting up base camp in the Euclidean 
world from which to tackle the hyperbolic world by 

• analytically extending the fundamental solutions of Laplace's equation, 

• computing the extended 5-sources, and 

• extracting Minkowskian propagators and sources from these extensions. 

Because we will need the fundamental solutions in M. 3 as well as M 4 , we work in this section 
with M. n for n > 3. The fundamental solution |T96j G n for Laplace's equation in M. n is 

I r 2-n 

A n G n (x) = -8 n (x), G n {x) = -, n > 3, (50) 

uo n n — I 

where uj n = 2n n l 2 /T(n/2) is the area of the unit sphere in W 1 and 

r(x) = Vx 2 

7 I was unaware of this work when developing acoustic and electromagnetic wavelets, and thank Lou Felsen 
for pointing it out. 
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is the Euclidean distance. 8 To extend G n analytically, we need only extend r. Define the 
complex distance from iy (source point) to x (observation point) as 

f(x — iy) = \J (x — iy) 2 = a/ r 2 — a 2 — 2ix ■ y. (51) 

Fixing y 7^ 0, the branch points form an (n — 2)-sphere in the hyperplane y ± of R n , 

S(y) = {x Gf" : r = a, x ■ y = 0}. (52) 

We will use the cylindrical coordinates 9 (p, cr,£) given by 

£ = yx, p= yjr 2 -^ (y = y/a), (53) 

with <r a vector in the unit sphere S(y) of the hyperplane y ± . The intersection of S(y) with 
the half-plane through y and cr is the point (a, cr, 0). If a; follows a simple loop surrounding 
this point (i.e., if the observer links the circle), then f changes sign. To make f single- 
valued, we must prevent the completion of such loops by choosing a branch cut consisting 
of a hypersurface with S as its boundary. The branch cut must be chosen so that f reduces 
to the usual distance on W 1 , 

y — > =>- f{x — iy) — > +r(cc). 
The simplest such cut is obtained by requiring 

Re f > 0, 

which gives the disk spanning S(y), 

V(y) = {x : r < a, x ■ y = 0}, &D = S. 

The most general branch cut is a membrane obtained from T> by a continuous deformation 
leaving its boundary intact. 

Note: By 'branch cut' we really mean a slice of the branch cut of r(z) at constant y, since 
the source is taken as fixed. 

Fixing y ^ 0, write 

f = yjx — iy) 2 = p(x) — iq(x) (54) 

with the t/-dependence understood. Then p(x) > in R n , T> is characterized by p = 0, and 
5 by p = q = 0. Squaring (jo^j) gives 

r 2 — a 2 — 2ia! ■ y = p 2 — q 2 — 2ipq 

2 2 2 2 / rr \ 

=^>r — a =p — g , x ■ y = pq. (55) 

8 The potential in K 2 is G-2 = —(2tt)^ 1 lnr. Although our method works here as well, this case is somewhat 
special and will be treated elsewhere |K0xj . 

9 To help visualize the geometry, the reader may think of the case n = 3, where S(y) is the circle forming 
the equator of the sphere whose north pole is y and the unit vector cr is represented by the azimuthal angle 
<j) in the plane y ± . 
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The cylindrical coordinates are thus given by 

of = pq 

2 2 2 2 2 2 

a p = a r — p q 

2/2 2 i 2\ 2 2 

= a (p — q + a ) — p q 

= (p 2 + a 2 )(a 2 -q 2 ). (56) 

In particular, note that \q\ < a — \y\, i.e., the imaginary part of r(z) is bounded by the 
modulus of the imaginary part of z. 




Figure 1. Plots of Re r (left) and Im r (right) with n = 2 and y = (0, 1). D(y) is now the interval 
[— 1, 1] along the x\ axis. The graph of Re f(x) is a pinched cone, and the jump of Im f(x) across 
D(y) is 2yT — xf, hence the cut is circular. Note from (|58[) that far from the disk, Re f ~ r and 
Im f ~ — cos 9 (since a = 1 here). 

It follows immediately from (|5fij) that the level surfaces of p(x) and q(x) are 



f P 2 ^ 2 1 

E„ = {constant p > 0} = < a? : ^ + — = 1 > 

I p z + cr p A J 

i?* = {constant < ±q < a} = {x : 2 P 2 - ^ = 1, ±f > J> . 



(57) 



a* — q* q" 



The Ep's are a family of oblate spheroids filling the complement of T> in M. n , with S as 
their common focal set, and the H^s are the orthogonal family of upper and lower semi- 
hyperboloids, also iS-confocal and joining in T>. As p — > 0, E p converges to a double cover 
of V, a fact that will be important in our computations. Similarly, as q — > 0, converge 
to the upper and lower covers of the complement of T> in the hyperplane £ = 0. Finally, as 
q — > ±a, the semi-hyperboloids collapse to half-lines: 

Hf = {q = ±a} = {±Xy : A > 0}. 

For the pulsed beams, the E p 's (with n = 3) will be wave fronts, the H^'s give the 
orthogonal surfaces of radiation Bow, and will be the forward and backward beam 
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axes. Use will also be made of the far zone approximation, where the observer is far from 
the disk. By (IHH1) . 

r ^> a =^> p pa r, q pa a cos 9 (58) 

so that E p becomes the sphere r = p and become the cones cos# = ±q/a. 

The complex distance thus provides a natural set of coordinates in R n , called oblate spheroidal 
(OS) coordinates, given by 

(p,q,cr) : p>0, -a < q < a, a E S(y). 



We now define the point source at iy by 





1 z2—n 


5 n (x - iy) = -A n G n (x - iy), 


Gn(z) = 


u n n — 2 



where A n is the (distributional) Laplacian in x. It can be shown FvOOj that for any y, 
5 n (x — iy) is a compactly supported distribution in x e M. n in the sense of Schwartz [T96 . 
Although the proof is somewhat involved, the supports are easily found. For even n > 4, 
G n (z) is analytic wherever z 2 ^ 0, in which case A n G n (z) = 0. Hence 10 

supp x S n (x — iy) = S(y) for even n > 4. (60) 

But for odd n, G n inherits a branch cut from r, where differentiating across the discontinuity 
contributes to the support of 5 n . Thus 

supp x 5 n {x — iy) = V{y) for odd n > 3. (61) 

(The same holds for n = 2 since 2irG2{z) = — lnf has a branch cut on the intreval T>.) The 
distribution S^(z) will be computed later along with its time-dependent version for pulsed 
beams. 

To illustrate the above, we work out the case n = 1 which, although trivial, is indicative as 
will be seen. Recalling our notation x for the sign of x, we have for z = x — iy G C 

f = \J (x — iy) 2 = x(x — iy) = \x\ — ixy, i^O. 

(Note that this is not simply the distance in C between iy and x, which would be y|x — iy\ 2 
and not yjjx — iy) 2 .) Since the 'unit sphere' in 1R consists of x = ±1, its 'area' is LO\ = 2, 
and ()50|) gives the correct solution in M: 

= - y => d 2 x d = -<Si(s). 

10 In the Minkowski space R"~ ' , y becomes the time axis, S a slice of the light cone r — t, T> a slice of 
the future cone r < t, and (|60|l and translate to Huygens' principle in even spacetime dimensions and 
lack thereof in odd spacetime dimensions; see |K00| . 
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Therefore 



G x {x - iy) = - X ^ , d x Gi(x - iy) = -~ + iy6i(x), 

hence 

~8 x {x - iy) = -d 2 x G x (x - iy) = 5 1 (x) - iy5[(x) (62) 



with Fourier transform 

~5 x {k,y) = \ dxe- ikx S 1 (x-iy) = l + ky. (63) 



4 Point sources in complex spacetime 

Following our 'Euclidean strategy,' we bring time into the picture by complexifying Euclidean 
spacetime: 

x E = (x, -u) e R\ y E = (y, t) e M 4 
z = x E — iy E — (x — iy, —u + it) E C 4 . 

This can be rewritten as a complex Minkowski vector 

z = x — iy = (z, ir), z = x — iy, r = t — iu, z 2 = z 2 — r 2 , 

where 

x = {x,it), y = (y,iu) G R 3 ' 1 
are regarded as real Minkowski vectors with pseudonorms 

2 2 j2 2 2 2 

x —x — t , y = y — u . 
Now consider the fundamental solution for the Laplacian in M 4 and its holomorphic extension, 

A 4 G 4 (x, u) = (A,,. + dl)G 4 (x, u) = -5 4 (x, u). (64) 

We have already seen G(z) = —G 4 (z) in (|46jl. but there it had the above form only in T + , 
being antisymmetric because of the sign factor u in fJ42jl . This 'twist' was the result of 
applying the AST to a spacetime field and is another example of the pitfalls of starting from 
the Lorentzian world. To be precise, u gives the orientation of T relative to the way M 3 ' 1 is 
embedded in C 4 (recall the one-dimensional case (|25)l). thus it has no role in the Euclidean 
world. In fact, we can now see that the restrictions to T ± of the extended relativistic field 
*&(z) in ()46|) are part of a single holomorphic function, just as uG(z) are both part of G 4 (z). 
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Before defining point sources in complex spacetime, or complex event sources, we must 
understand how real point sources in R 3,1 fit into this picture. Naively, it might be hoped 
that the substitution u — > it into gives a propagator for the wave equation, 

□G 4 (x, t) = (A;,, - <9 t 2 ) G? 4 (x, t) " - 5 3il (aj, t). (65) 

We will see that this fails for the following reasons: 

• G4 is singular on the light cone x 2 = and must be defined as a distribution in R 3 ' 1 by a 
limiting process. 

• When properly defined in M 3 ' 1 , it turns out to be sourceless. 

• Propagators are related to causality, which depends on the existence of the light cone and 
hence makes sense in M 3 ' 1 but not in M 4 . Thus G^{z), which comes from M 4 , cannot itself 
be a propagator. 

To find the extended propagators, note that 

z 2 = z 2 -r 2 = f 2 -r 2 = (f-r)(r + r), (66) 
which gives the partial-fractions decomposition already encountered in 



Ga(z) = iD~(z) — iD + (z), D ± (z) 



iiir 2 r(T =F f) 



(67) 



Since formal differentiation gives 

DD ± (z) = 0, 

D ± can have sources only at x 6 D or f = ±r. We will show that D ± {z) are the proper 
extensions of the advanced and retarded propagators in M 3,1 . Note that choosing the 'un- 
physical' branch of f simply reverses the sense of causality: 

f-v-r => D ± (z) -> -£> T (z). (68) 

Thus it will suffice to study the properties of _D + (z). Its time behavior is governed by the 
retarded Cauchy kernel: 

C(r-f) 



2m{r — f) 2tt (u — q) + i(t — p) ' 

which gives the time-domain radiation pattern |HY99I . This shows that the ellipsoids E p are 
wave fronts, ie., surfaces of constant retardation for D + . An observer fixed at x will see a 
pulse peaking at time t(x) = p(x), with duration 

T(x) = \u-q(x)\. 

The peak magnitude 
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is constant along the hyperboloids Thus, apart from the attenuation factor 1/r, the peak 
value of |.D + (;2)| remains constant along the hyperboloids. This shows that the radiation 
Hows along H+ if < q < a and H~ if —a < q < 0. In the far zone, C{r — f) has peak 
magnitude 

n(9) = \C(r-r)\ t=r = ~ -, r»a, (69) 

Z7r \u — acosu\ 

showing that there are three qualitatively different cases, depending on the causal character 
of the imaginary source point. 

1. If y is timelike (y 2 < 0), then \u\ > a and D + is a smooth pulse outside the source 
region T>. Furthermore, the peak radiation pattern 7Z(6) is an ellipse with eccentricity 
a/\u\ and the source at one of the foci. If y is in the future cone (u > a), then the 
semimajor axis of the ellipse points in the direction y, so that the exploding wave 
D + (z) is emitted along y. It will be shown that G 4 = iD ~ — iD + is sourceless, so 
D~(z) is an imploding wave that is simultaneously absorbed along —y. 

If y is in the past cone (u < —a), then D~(z) is absorbed along y while D + (z) is 
emitted along — y. As y 2 —* 0~, the ellipses become more and more eccentric and the 
pulsed beams become sharper and sharper. Notice that they have no sidelobes, 11 hence 
could be useful in applications such as radar tracking and directed communications 

[K96llK97llKnT| . 

2. If y is future lightlike (u = a), then D + is singular on the ray along y at t = r. If y 
is past lightlike (u = —a), D + is singular along — y at t = r. In either case, the peak 
radiation pattern 7Z(6) is parabolic. 

3. If y is spacelike (y 2 > 0), then \u\ < a and D + is singular at t = r on the cone 
cos# = u/a, with a hyperbolic radiation pattern. 



Only Case 1 gives reasonable pulsed beams with a chance to be realized with finite energy, 
although Cases 2 and 3 should also be of interest since G±(z) is holomorphic for all z 2 ^ 
0. (Quantum field theory comes to mind, where virtual energy exchange can occur along 
spacelike directions.) We therefore assume from now on that y G V±, so that z belongs to 
the causal tube, 

zeT = T + UT_, (70) 

as already introduce in (|13J) starting from Minkowski space. 

11 Sidelobes are angular patterns resulting from interference between waves coming from opposite ends of 
the source. They can cause problems in radar and communications. 
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Figure 2. Time-lapse plots of \D + (x — iy)\ in the far zone, showing the evolution of a single pulse 
with propagation vector y = (0,0, l,iu). We have taken the slice X2 = 0, so that the source disk 
becomes the interval [—1,1] on the xi-axis and the pulse propgates in the X3 direction of the X1-X3 
plane. Clockwise from upper left: u = 1.5, 1.1, 1.01, 1.001. As u — > 1, y approaches the light cone 
and the pulsed beams become more and more focused. 

We are now ready to learn how D ± (z) are related to the wave propagators in M 3,1 . It suffices 
to focus on D + (z). Fix y £ V + and define the partial boundary values 

D + (x±i0) = \imD + (x±iey). (71) 

The limits in (171)) do not depend on the choice of y £ V + , hence we may choose y = (0,iu) 
with u > so that 

D + {z) =D + (x,r) = — — , r = t-iu, 

St7T z r r — r 

and the Plemelj jump conditions [T96J give the distributional relations 

1 1 1 1 S(t-r) 

8in 2 r e\o t — r =p ieu Sm I r t — r 8wr 

where V is the Cauchy principal value. If we require Huy gens' principle jBC87j to be valid in 
M 3 ' 1 , then the principal value terms must be eliminated and the only combination acceptable 
as retarded propagator is 

5(t - r) 

D + ( y x) = D + ( y x~i0)-D + (x + i0)= \ ■ ( 72 ) 
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This combination of boundary values, already encountered earlier, will be called the Minkows- 
kian limit. It does indeed give 



DD + (x) = -6 3 (x)5(t) = -5 3il {x). 



(73) 





t igure 3. \D + {x - iy)\ 2 with y = (0,0,1,1.01*) in the near zone immediately after launch, 
evolving in the x%-x 3 plane with x^ = as in Figure 2. Clockwise from upper left: t = 0.1, 1, 2, 3. 
The ellipsoidal wave fronts and hyperbolic flow lines are clearly visible. The top of the peak is cut 
off to show the behavior near the base. The spikes represent the branch circle, whose slice with 
X2 = consists of the points (±1, 0, 0). 
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The substitution r — > — r as in (jb8j) now gives 

6{t + r) 



D-{x) = D-(x-iO) - D-(x + i0) 



Attt 

OD-{x) = -5 3>1 {x). 

Equation (fTHJ) is the desired hyperbolic counterpart of Its derivation confirms the 

points made about the failed attempt 

• The distributional limit (ff2^) played a key role in reproducing Huygens' principle. 

• By (lf)7j) . the Minkowskian limit of iG±(z) is 



R(x) = iG A {x-iO) -iG^x + iO) = D + (x) - D~(x) 
_ 8(t - r) 5(t + r) 
47ir Airr 



(74) 

This is the Riemann function |T96j . which solves the following initial- value problem of the 
sourceless wave equation: 

DR(x) = 0, R(x, 0) = 0, d t R(x, 0) = 5 3 (x). (75) 

iG±(z) can therefore be considered the extended Riemann function. It is also identical with 
Synge's elementary wave function [S65J. 

• Causality has no meaning for G±(z) and appears only when a branch cut is chosen for f, 
as Equations (JHZj) and (JtjBj) confirm. 

The limits fl?2^) and (|74^) are typical of hyperfunction theory [K88l II92j , where distributions 
are represented as differences of boundary values of functions holomorphic in 'local' wedge- 
like domains surrounding the support. In general, there is no preferred set of such domains 
and it is necessary to use sheaf cohomology, which makes the theory rather abstract. In our 
case, however, the two domains T ± suffice due to the natural cone structure of relativistic 
equations. 

Following (f73|) . we now define the point source at iy as 



S 3 ,i(x -iy) = -D x D ± (x - iy), y 2 < 



(76) 



We list some of its basic properties. 

• The left side of f!76|) is independent of the sign on the right. Recall that \3G±(z) = 
wherever G4 is holomorphic. For z = x — iy G T, 

z 2 = =3- x 2 = y 2 < and x ■ y — 0, 

but x ■ y cannot vanish since both vectors are timelike. Hence G±(z) is holomorphic in T and 

DD + (z) -DD-(z) = DiG 4 (z) = 0. 
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• The Minkowskian limit of S 3j i(z) is 5 3} i(x). By definition, 

£3,1 (x) = S 3j i(x - z'O) - S 3)1 (x + z'O) 

= -DD + (x - z'O) + QD + (x + z'O) 

= -D J D + (x) = -5 3 ,i(x). (77) 

• ^3,1 ( x ~~ w) JS supported in the world tube swept out in M 3 ' 1 by V(y) at rest: 

supp 6 3 ,i(x - iy) = {(x, it) : x e V(y)} = V(y), Vy e V± . 
This follows since D-D + = outside the singularities and 

- Im (r - f) = u - q Vz G T, 
so the only singularities come from the 1/f factor. 

Note that while the Minkowskian limits D ± (x) and £3,1 (a?) are Lorentz- invariant, the ex- 
tended propagators D ± (z) and their sources S 3) i(z) are frame-dependent, the preferred frame 
being the rest frame of T>. This can be traced back to the fact that we have obtained the 
factorization and associated splitting (|57jl by choosing a branch cut for f(z) in a par- 
ticular Lorentz frame. Of course, since G±(z) is Lorentz- invariant, we may choose to do the 
splitting in any other frame. 

5 Driven complex sources 

Suppose we 'drive' a point source fixed at iy with a real time signal go(t). The resulting 
retarded wave is the convolution 

00 



W(z,r) = I dt' D + (z,r-t')g (t 

-00 

1 f°° g (t')dt' 



\m 2 f ./-oo t — V — f 



( ' (T r) =-g(r-r)G 3 (z), (78) 



where G 3 (z) is the holomorphic Coulomb potential |N73L iKOlaj and 

is the analytic-signal transform (|21|) of g . The source distribution for the associated field 

W(z,r) = 9 -^1 (80) 

is defined by 

S(z) = -D x W(z). (81) 
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Since g is analytic off the real axis and 

— Im (r — f) = u — q ^ Vz e T, 

it follows that W(;z), like D ± (z), is analytic outside the world tube T> swept out by T>(y) at 
rest. Moreover, formal differentiation gives 

DW(z) = 0, 
therefore S(x — iy) is also supported in x e T>. 

Examples: We give three driving signals that will be needed later, with their ASTs and 
radiated waves (jSUjl: 

g {t) = S(t) 9 ^ = ^ W = D + (z) (82) 

9o(t) = 1 g{t- iu) = u/2 W = -(u/2) G 3 (z) (83) 

9o(t) = e-™ g(r) = 6{u, u) W(z) = g(r) B„{z), (84) 



where C is the Fourier transform of the Cauchy kernel 

C(uj,u) = u&(uju)e- UJU (85) 

and 

B ^ = (86) 

is the time-harmonic complex- source beam, applied widely in engineering |HF01j . In the last 
equality of we use d 

\q\ < a =3- Sgn (u — q) = Sgn u =3- Q(u(u — q)) — Q(uu). (87) 

Note that to make (J83|) a special case of (J82)) with uj = 0, we must define 0(0) = 1/2. 

The far-zone approximation (|58|) to (|84j) shows that W is indeed a beam in the direction of 
u y which becomes more and more focused as y approaches the light cone: 



r>a W{z)^ue(uou)e- u){u - acosl 



e -iw(t-r) 

Airr 



This is not surprising, since (|H4"|) is a Fourier component of the pulsed beam (I82J) . The beam 
W is exponentially stronger for q = u than for q = —u, so the beams are directed. Equation 
(JHH|) shows that the imaginary retardation u — >• u — q serves to focus the beam. 

Just as the spatial displacement of a point source from x = to iy expands it to a disk, 
the temporal displacement of the impulse from t = to iu gives it duration, as seen in 
). We may interpret the parameter \u\ as a response time for the source at iy. A large 
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response time suppresses rapid variations in the driving signal. This is modeled 12 by the 
factor Q(uu) e~ u ' u , which acts as a filter to smooth the signal. The larger the source, the 
longer the response time since the excitation can travel to different parts of the source before 
emitting a wave. Thus, a rough but intuitive way of understanding the timelike character of 
y is to note that the time needed for a signal to travel from the center to the rim of the disk 
is a, and therefore \u\ > a. 

With the interpretation of u as a 'response time,' all four components of y have a direct 
significance in terms of the source itself, without reference to the radiated beam. 



6 Main results on sources 

In this section we present our main results on the source distribution for beams W(z) of the 
general type (JHUJl . These will be proved rigorously in the Appendix. The choices ()82|) -([5I j) 
for g yield the sources for pulsed, static, and time-harmonic beams generated by complex 
event sources. Furthermore, choosing a plane wave for the test function will yield the space- 
time Fourier transforms of the sources and their beams, giving valuable insight into their 
propagation properties and making them a potentially useful computational tool. 

For clarity, the results are stated as theorems. However, I have included some discussion 
to help make them digestible for readers without training in the art of arid mathematical 
discourse. A more extensive discussion of the results and their interpretation is given in the 
next section. 

To compute S(z) as a distribution, we must deal with its singularities. This will be done 
by shielding T> with an ellipsoid E £a and taking the limit e — * once the computations are 
complete. Thus let e > and 



W £ {z) = Q{p-ea)W{z) 




if x is outside E t 
if x is inside E, 



eai 



The singularities on T>, consisting of the discontinuity in the interior of the disk and the diver- 
gence on the boundary S = &D, have been replaced by a uniformly finite jump discontinuity 
across E ea . The regularized source, defined by 



S £ (z) = -DW e (z), 



is therefore supported on the world tube swept out by E ea at rest, 

E £a = {x = (x,it) : x G E £a }. 
The 'bare' source S will then be defined by 



S(z) = 1imS e (z), (89) 

e\0 



12 Admittedly this model leaves something to be desired since g decays slowly even if go has compact 
support. In particular, it is not even causal: the vanishing of go(t) for t < does not imply the same for 
g(t — iu); in fact, the latter cannot vanish on any interval because it is analytic. One way to improve the 
model is to represent the response by derivatives of g(r), as in (|48J) . This suppresses the low frequencies and 
amounts to putting go through a 'band-pass filter' of the form cj Q 0(a;u)e~" u centered around u) = a/u. 
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where the limit is taken in the distributional (weak) sense by 'smearing' over test functions. 

Note: The definition ()89|) is necessary on conceptual as well as technical grounds for the 
following subtle reason. We have chosen a branch of the complex distance f to reduce to the 
usual positive distance as y — > 0, but nowhere have we actually enforced this in our equations 
- until now! The definition ()89)) clearly communicates our choice to the equations. 

Theorem 1 (Shielded complex source) 

For given timelike y, S £ (x — iy) is a (Schwartz) distribution [T96J in x E M 3 ' 1 supported on 
E ea . It is regular in t (no smearing necessary) and acts on a spatial test function f(x) as 
follows: 



SeJ) 



a a 
~2a~ 



gf 

if 



a*a f a dq ° 
+ — / — gff 



p=ea 



(90) 



where a = ea — ia and a* are the north and south poles of E ea and we used the notation 

f27T 



/(f) = f(p,q) = — J f(p,q, 

2 f- = 2drf = dj + tdj 

without implying analyticity in f = p — iq. S e is given directly (without smearing) as 



(91) 



S e (z) = W{z) { i5{f -a)- i5{? - a*)- 



a 
r 



5(p — ea) (dp + id q ) 



(92) 



where 



8(f — a) — 5(p — ea)5(q — a), S(r — a*) = 5(p — ea)5(q + a) 

are point sources at the north and south poles of E ea . The normal derivative d p is interpreted 
as a double layer, and the tangential derivative d q as a surface flow on E Ea . The entire 
distribution is modulated in space and time by the values of the beam W(z) on E £a , including 
the retarded analytic signal g(r — r) of its driving function responsible for collimating the 
beam. The distribution (192)1 must be applied as a differential operator to a test function, then 
integrated. 



Theorem 2 (Unshielded complex source) 

For given timelike y, S(x — iy) is a distribution in x e M 3 ' 1 , regular in t, acting on a spatial 
test function f(x) by 



S J) = g(r,a)f(0) +ia — g(r,q) f p {iq)+if q {iq) 

Jo Q 1 



(93) 
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where 



(94) 



is t/ie average of g(r — f) over t/ie jitmp across D . The integral converges because the conti- 
nuity of f(x) across D and its differentiability on T> imply that the partial derivatives of f 
on T> are 0(q): 



(95) 





fp(iq) = - h > f q ( i( i) = -- 

CL jD 


f P 




where f p = d p f, f% = d^f . In cylindrical coordinates, with 


q= a/o 2 - 


P 2 , 




(SJ) = ~g(T,a)f(0) + r^~g(r,q) 

Jo q 


af P + iph 




S can be expressed without smearing as 




S(z) = ~g(r, q) U(x) + ^^l^ (afy + iffy) 


I 



(96) 



(97) 



Therefore, S(z) has a real point source at the origin and a surface distribution on T> consist- 
ing of a radial flow d p and a double layer d^. The entire distribution is modulated in space 
and time by g(r, q). 



Remark 1: Constructions of complex sources exist in higher dimensions [KOx], but only for 
n < 3 does 5 n (z) contain a real point source. (See Equation (|52Jl for the trivial case n = 1.) 
The reason is that for n > 3, S inherits the factor p n ~ 3 from the volume element, and this 
vanishes at the poles of E £a . 

Remark 2: The point source 8s(x) in ()97j) is not apparent in Figure 3, which seems 
puzzling. However, recall that this point source is the limit of the two point sources in (|92j) 
and, as explained below (|185|) . these point sources are 'removable singularities' resulting 
from an integration by parts. Without them, S e (z) has the equivalent 'local' expression 
(|186j) . However, that expression does not have a limit as e — > and must be regularized by 
subtracting and adding a Taylor approximation to the numerator, as was done in |K00j for 
the general case of lR n . (This procedure gives a generalization of the Cauchy principal value 
to arbitrary order; see jGS(i The resulting limit distribution (S,f) is nonlocal, i.e., it 
cannot be written in unsmeared form like (J97|) . We write 'local' is in quotes because (J97|) still 
contains differential operators. However, these can be interpreted locally as double layers 
and flows, which would be impossible if ( S, f ) had subtractions. This is especially important 
if one is interested in the possibility of realizing such sources, i.e., building instruments that 
can emit and, by reciprocity, detect pulsed beams of the above type. 

Remark 3: An even more intriguing possibility is that some naturally occurring fundamen- 
tal processes can be associated with emission, absorption, or radiation by complex source 
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points. Natural extensions of the present construction exist for Maxwell's equations |K02j . 
and even for massive Klein-Gordon and Dirac fields |K9QI , although the sources in the latter 
cases have not yet been studied. In fact, one of the first applications of what I have called 
complex distance occurred in the derivation by Newman et al. of charged spinning black 
holes; see jN65llN73j . and the recent work [KOTallNOl] . 

Corollary 1 (Event sources, static and time-harmonic point sources) 

Letting go(t) = 5(t) as in (J82j) gives g(r) = 1/2hit and 



W(t, q) 



T 



2tt(t 2 + a? - p 2 



2lTZ 2 



z G V. 



Therefore the event source at iy is 



2lTZ 2 



2ixqp 



(9* 



Letting go(t) = 1 as in (}83|) gives the static point source at iy: 



Uz) = 8 3 (x) + 8( ° P)m (ad p + ipdi) 
27rq p 



(99) 



Letting go(t) = e tujt as in f!84|) gives the analytic signal for the time-harmonic point 
source at iy: 



g(r,q) = uQ(uu)e luJT cosh(a;g). 



(100) 



Much valuable information is gained by computing the Fourier transforms of the pulsed 
beams and their sources. Although this seems at first a difficult task because of the convo- 
luted spatial dependence of the beams through f and its oblate spheroidal coordinates, the 
result is intriguing and 'unreasonably' simple. 

Recall our notation for Minkoswki space M and its dual Fourier space M'. Given a spatial 
direction y, we use the cylindrical coordinates 

x = (p, £, it) k = (h, I, iuj), p ■ y = h ■ y = 

£ = y ■ x p = |p| = \/ r 2 — £ 2 r = \x\ 

l = yk h=\h\ = Vk 2 -1 2 k = \k\ (101) 

so that / and h are the longitudinal and transverse wave numbers with respect to y and the 
pairing between M' and M is 

k-x = k- x~ut = h- p + l^ — ujt = hpcoscf) + Z£ — uot. (102) 
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Theorem 3 (Fourier transform of shielded source) 

• Given y 6 V with \y\ = a > and e > 0, define the complex wave vector k £ 



uo F = lu — iel 



K = \v\ h 



7] = e — i = a/a 



L = I — ieuj 



which satisfies 



% = \ v \ 2 k 2 = k? + e z k 



2;„2 



(103) 



(104) 



and thus preserves the complex light cone. Then the Fourier transform (|18|) with respect 
to x of the shielded source S e (x — iy) is 



S e (k, y) = g(uj, u) e l£uja Q(k E , y), Q(k e , y) = cos(/i e a) + — sin(/i e a) 



(105) 



where 



and 



g(u, u) = u Q(uu) e uu g (uj) = C(u, u) g Q (u) 



(106) 



(107) 



The transformation k*—>-k £ is the product of a scaling k i— *■ \rj\k and a real rotation in the 
l-iu plane, or imaginary Lorentz transformation in the 1-uj plane, given by 



I — eiuj 
\flTe 1 



ILO I 



iuj + el 

7TW 



Q plays the role of a focusing Biter off the light cone, and on the light cone it simplifies 



as 



k 2 = Q <-> l = ±tfi <-> n(k,y) = e ±i ^ a . 



(10* 



Since \i is imaginary for propagating waves and real for evanescent waves, Q amplifies 'for- 
ward' waves and dampens 'backward' waves. It does not depend on the branch of the square 
root since it is even in fi. 



Corollary 2 (Fourier transforms of bare sources) 

The Fourier transform of the bare source S(x — iy) is 



S(k, y) = g(u, u) Q(k, y) = g(u, u) { cos(fia) H — sin(ita) 

/i 



(109) 
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where 



y. = Vh 2 - u 2 = Vk 2 - I 2 



can be real or imaginary. The Fourier transform of the event source at iy is 



$3,i(k, y) — C(u, u) Q(k, y), C{uj, u) = u Q(uiu)e 



(110) 



and that of the point source at iy is 



5s(k,y) = cos(ha) + — sm(ha). 



(Ill) 



Amazingly, complex sources are much simpler in Fourier space than in space or in spacetime! 
This shows that effective computations can be performed with them and their radiated 
beams using 'fast' numerical methods like the FFT. For example, spacetime convolutions of 
complex sources with arbitrary 'densities' can be performed with ease by multiplying the 
Fourier transforms. 

The surprising simplicity and symmetry of the above expressions may be a 'miracle,' but it 
cannot be an accident. Its origin and some consequences are discussed in the next section. 

Note that the expression (jlllj) is actually valid for 5 n (k, y) with any value of n > 1, as will 
be proved elsewhere. In particular, note that it holds in the trivial case n = 1 (|63|). where 
a = y, I = k, and h = 0. 

7 Interpretation and discussion of results 

Let us attempt to understand some of the expressions given in the last section. Our discussion 
is necessarily somewhat speculative, undertaken with the desire to add qualitative value to 
the raw mathematical equations. 

• The bare source S(z): Looking at Equation (|HTj) . 



is the retarded propagator for the wave equation in two space dimensions, with a = \y 
playing the role of time T!)(i : 




(112) 



we note that 




e(a-p) 



(A 2 - dl 



)G 2 , 



i — — 



5(x 1 )5(x2)S(a). 
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The origin of this analogy is not difficult to find. By complexifying the distance function 

r(x) — > r(x — iy) = a/ r 2 — a 2 — 2ix ■ y, 

we have implicitly introduced a as 'time' measured along the axis defined by the unit vector 
y, in exactly the same way as complexifying the Euclidean time u opened up the light cone 
and gave rise to the physical time t. In a coordinate system where y = (0, 0, 1), the branch 
cut T>(y) is simply a slice of the 'future cone' 



V+(y) = {{x u x 2 ,is) : p = \jx\ + x 2 2 < s} C M 2 ' 1 

at s = a, with p as the distance function in the orthogonal 'space' M 2 . The Heaviside function 
0(a — p) thus merely enforces 'causality' and the factor y/ 'a 2 — p 2 in the denominator reflects 
the absence of Huygens' principle in R 2,1 , where an 'observer' at (xi,x 2 ) will 'hear' nothing 
for s < p, then a 'sonic boom' at s = p with a decaying tail for s > p. 

Does this give any insight to the nature of S(z)7 It suggests looking at a as an evolution 
parameter. Think of 'morphing' from a point source to a disk source by gradually changing 
a, then ()112|) shows how the source flows while evolving. 



The shielded source <M 



S e (z) = W(z)U5(f-a)-i5(r-oc*)- - 5(p - ea) (d p + id q ) \ . (113) 



a 



r 



2 



This expression is extremely simple, reflecting the regularity achieved by replacing the singu- 
lar disk D with the oblate spheroid E ea . The first two terms are a pair of real point sources 
at the north and south poles of E £a , 

a = ea — ia = rja, a* = n*a. 

As already mentioned, the terms with the normal derivative d p and the tangential deriva- 
tive d q may be interpreted as a double layer and a flow on E Ea , although this needs some 
discussion. Further insight is gained from the Fourier transform ()105|) . 

SJk, y) = g(u, u)e iea Vt{k £ , y), Vt(k, y) = cos(pa) + - sm(pa). (114) 

p 

In the proof of (jll4j) in the Appendix, we saved the final details for this discussion because 
they shed light on the nature of the source, and also because they contain some spectacular 
cancellations and reveal an amazing hidden structure, namely the complex mapping k — » k £ 
of Fourier space associated with the replacement of the disk source D by the Huygens source 
S £ that generates the identical field outside E £a by emitting 'secondary wavelets.' 

By (HHHD, (HEHD and (JMI, 

Q £ (k, y) = I -h + ieh + V*vh 
Io = cosh(u; e a) — ie sinh(a; £ a) 

l\ = cosh(cj e a) — cos(p £ a) (115) 
I 2 = sinh(^ e a) - [yj £ j p £ ) sin(/i e a) 

h = (l/Ve) sm(p £ a). (116) 
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The derivation shows that Jo is due to the point sources at a and a* and the other terms 
are due to the the combined double layer and flow. Specifically, Ii and I 2 come from the 
p-derivative and I3 comes from the ^-derivative. The first terms of Jj and I 2 entirely cancel 
the point-source term I , resulting in 

tt e {k, y) = cos(jU £ a) - ie(u £ /fi £ ) sin(/i e a) + (e 2 + l)(l/fi £ ) sin(/x e a). 

This cancellation is related to the fact that the point sources in J resulted from an integration 
by parts, which in turn was needed to get a local expression for S(z) (see the note below 
IflBSjl ). But 

(e 2 + 1)1 — ieuj e = (e 2 + 1)1 — ie(uj — iel) = I — ieuj = l £ 
which gives the final form ([1140 . 

The cancellations and simplifications taking place to yield this simple result appear to be 
'miraculous.' This could be merely good fortune or, more likely, an indication that the Fourier 
sources and possibly also the unexpected complex mapping k — > k £ , are more 'fundamental' 
than the spacetime beams we started with and should therefore be thoroughly understood. 

On the practical side, simplicity in the Fourier domain usually means enhanced analytical 
power and the existence of efficient implementations by 'fast' algorithms. The above Fourier 
sources offer a promising new tool, modeling processes of directed emission and absorption 
that could be useful in quantum as well as classical physics. 

As a first step, and an example of working with the Fourier sources, we compute the pulsed 
beam W(z) from its source. 



8 Weyl representations for complex-source beams 

How does the simple Fourier source S(k, y) radiate a pulsed beam as complex (!) as W(z)l To 
learn the answer, we now compute W from S(k, y). In the process, we prove a generalization 
to beams of Hermann Weyl's representation of time-harmonic spherical waves. 

From Equations (jHUj) . ([STjl and ()109|h recall our path from W to S: 

W(Z,T)= J < T - ? > 



Airf 

S{x-iy) = -U x W{x-iy) (117) 



S(k,y)= / dx e %k ' x S(x - iy) = g{uo,u)VL(k,y) 
Jm 

Q(k, y) = cos(/xa) H — sm(fia) 

where g is the AST of a driving signal go(t) obtained by convolution (j2Tj) with the Cauchy 
kernel C(r), so that 

g{oj, u) = C(u>, u) goiyj), C(u, u) = uQ(uu) e _um . 
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Formally, ()117j) implies 



W(x — iy) 



ctk e 



M> 



k-x y) 

k 2 



:ii8i 



but the right side must be defined since k 2 vanishes on the light cone. In spacetime terms, 
this corresponds to the fact that "initial values" must be specified in order to solve (jll7|) 
for W . More precisely, since we are dealing with — oo < t < oo, we need the behavior of W 
as t — > ±00, which describes the causal relation between the source and the solution. This 
amounts to a choice of contour in Fourier space that avoids the light cone singularities, so 
that the solution can be computed by residues. 

The temporal Fourier transform of W is 

dt e iuJt g(t - iu - f) = g(u,u)B u (z), 



Airr 



where 



B»(z) = — , (119) 

is the time-harmonic complex source beam (|86|h and translating the integration contour by 
t —>■ t + f, or t —>■ t — iq, is justified because g is analytic off the real axis and \q\ < \a\ < \u\. 
Therefore, we need to establish that the function 



U(z,u) 



dk e 



k - x V{k, y) 
k 2 



hdhdchdl e ihpcos *e il t ■ 



fX 2 + l 2 



hdh J (hp) I dl e a Z ■ 



(120) 



is identical with B u (z), provided the integration contour is chosen to give the known behavior 
of W. Note that 

Q(k, y) = — {// cos(//a) + I sm(fia)} 



2/i 



{(/i-z/)e^ a + (fi + il)e- ilMa } 



giving 



k 2 



+ 



2/x(/x + il) 2/x(/x — il) 



:i2v 



This key identity will reveal how the 'focusing filter' Q amplifies forward waves and suppresses 
backward waves. Inserted into ()120|) . it gives 



U(z,u) = 



hdh 



J (M / dl e m 



-ifia 



+ 



jj, + il pi — il 



(122) 



35 



The choice of contour thus amounts to picking a branch of /x = Vh 2 — uj 2 . Since the denom- 
inators \i ± il will give residues at I = ±z/i, the plane waves e l ' will be propagating when 
/i 2 < u; 2 and evanescent when /i 2 > a; 2 . Define the branch 



-is/up- 



h 2 



Vh 2 



h 2 <uj 2 
h 2 > u 2 . 



and note that all its values can be shifted to the right half-plane by the infinitesimal transla- 
tion \i — > n + 0. If we use this branch in (jl22j) . then by closing the integration contour in the 
upper or lower complex half-plane, depending on the behavior of e l1 ^, and using Cauchy's 
theorem, we obtain 



-i/ia 



di 



-oo 
oo 



(Tl 



+ 



e(f) e v« e -* 



6(Oe 



( = £ — ia 



e(-e)e- 



[A — il 

This gives U as the even part of a function U + defined by 



U(z,u) 
U + (z,u) 



U + {z,u) + U- 
roc hdh 



[p,0 



0(0 



2/i 



Mhp) e-^. 



(123) 



Thus, using the branch fi results in all waves propagating in the +y direction > 0) being 



amplified by e 



tpa 



\/w 2 — h 2 c 



and all those propagating in the — y direction (£ < 0) being 



suppressed by its reciprocal. The evanescent waves decay in both directions, as they should. 
Note that U (z, to) is even in z, but not in x alone. That is, while it has a preferred direction 
in real space, it does not have one in complex space. This is obvious since a pulsed beam in 
the — y direction will have the same behavior in —x as the original one has in x. 

Now recall from the discussion below (|69J) that the retarded pulsed beams propagate in 
the direction of uy, i.e., along y if u > and along — y if u < 0. But the factor Q(uu) 
in the Fourier transform of the Cauchy kernel forces the signs of u and to to be identical 
on the support of W. Therefore we need the growing exponential e* Ma associated with the 
+£ direction if uj > and the — £ direction if u < 0, and the decaying exponential e _va 
associated with the — £ direction if u > and the +£ direction if uj < 0. 

This shows that (jl23j) gives the correct value for uj > but the wrong one for u < 0. The 
correct branch of /x for ui < can be obtained by noting that r(z*) = r(z)*, hence 



Bjz*) 



Aitr 



BJz). 



Therefore we define U for negative frequencies by 



U(z,-u) = U{z*,u)* = U-(z,lu) + U-(-z,uj) 

"°° hdh 



U-(z,iu) = U + (z*,ujy = Q(0 







2/i* 



Mhp) e~^, 
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where p* is the branch 



Woo 2 - h 2 , h 2 < uo 2 
Vh 2 - oo 2 , h 2 > uo 2 . 



p = 

With this, proving ()12()|) reduces to the following. 



Theorem 4 (Generalized Weyl formula) The time- harmonic complex- source beam B u (z) 
has the following angular spectrum representation: 



Aur 



'U + (z,oo), oo > 0, £ > (large right component) 

U + (— z,uo), oo > 0, £ < (small left component) 

U + (z*,oo)*, uo < 0, £ > (small right component) 

k U + (— z*, uo)*, oo < 0, £ < (large left component) 



(124) 



where, for z = (p, () = (p,l; — ia), oo > 0, and £ > 0, 

U + (z,oo)=[ ^e* h -^ = U p :jz,oo) + U:jz,oo) } (125) 
with the propagating and evanescent parts given by 

^^[^h^^ (126) 
u ^=[^w J ^^- (127) 

The components U + (z,oo) and U + (—z,oo) in (|124j) are analytic continuations of one another 
across the plane £ = 0, with equal boundary values on that plane for p > a. The jump 
discontinuity, due to the branch cut T>, is imaginary and given by 

J(p, oo) = lim {Bu(p, e-ia) - B w (p, -e - ia)} 

e\0 

Q(a — p) , , r-r 

= % o HT—2 ^sh(oo^a 2 -p 2 ), (128) 
27ry or — p 

with spectral decomposition 

J(p,oo) = i / dhe*'*^ = i rhthMhp) 8 ^-. (129) 
Jr 2 ^ Jo ^ 

Remarks 

• In the limit a — > 0, ()124|) becomes Weyl's angular-spectrum decomposition of the funda- 
mental solution of Helmholtz's equation (see MW95 , pp. 120-125, where m = if/,), 



Airr 



Jr 2 ^fi 
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• Since f depends only on the squares of the components of z, it suffices to reverse only the 
component £ = £ — ia in the y direction. However, the beams are cylindrically symmetric 
about the y axis, therefore 

B u (p, -C) = B u (-p, -C) = B w (-z) 
and reversing ( is equivalent to reversing z. 
Proof. 

Equation ()124j) is a direct consequence of Formula (26) on page 9 of |E54j . With 

x = h, b = u, y = p, a = ( = (,— ia, £ > 
VV + a 2 = vV + C 2 = f, 

it states that the function 

'±i^h/(u 2 - h 2 ) c^W-v , 0</i<a- 



has Hankel transform 



^h/(h 2 -u 2 )e~^ k ' 2 -"\ u<h<oo 



dh f(h)J Q (hp)y/hp=^- , (130) 

r 



which gives ()124j) for positive uj and £. The other cases follow from the symmetries discussed 
above. 

That U + (p, — £ + ia, uj) is the analytic continuation of U + (p, £ — ia, a;) to £ < follows from 
the known analyticity of B^{z) outside of the branch cut T>. Recall that 

c ~ / 2 i (c ^2 /=FV a2 -P 2 ' °<P< a 

[VP -a , P>«, 

hence 

±ie ± "^ a *- f * /Any/a 2 - p 2 , 0<p<a 



'yf? 11 * finyfi? - a 2 , p>a 



B^{p,i-ia) 

and the jump across £ = is indeed given by ([128)1 . By ()125|) . the boundary values of the 
propagating and evanescent parts are 

hdh 



\imU + (p, e — ia,uj) = i I /"""" JJhp) e aVuj ' z - h ' 2 = iA 
lim U p :jp, ia-e,u,) = i J (hp) e -^^ = iA ' 

e\0 Jo 2V^ — fi 

lim tO, e - ia, uj) = f° ° J^L= J (hp) e™^? = B + iC 



lim[/+Jp,m-£,u;) = / J (fep) e -»°v/^ = B - iC 
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with A, A', B and C real. Hence the real part of B w is continuous and the jump across £ = 
is 

J{p,uo) = iA- iA' + 2iC 

= i [ - === J (M sinh(aV^ 2 - ^ 2 ) + 2 / 7== = Jo(hp) sin(aVh? - u 2 ) 
Jo \Jijj l - h 2 J w \Jh l - u 2 

Jo(hp) sm(pa). ■ 



Now that we have the correct contour, we also know the Fourier transform of the pulsed 
beam. By (fT2Tjl . 

p z + l z 2p [p + O + il p + — il) 

9 Electromagnetic wavelets revisited 

The scalar complex point sources will now be used to construct sources for electromagnetic 
wavelets. The most direct formulation is in terms of Hertz potentials, which are reviewed 
briefly. 

Although Hertz potentials have a long history in electrodynamics, they have been sadly 
ignored in most modern textbooks. Many books that do mention them have only a short 
section on the subject, usually under various specialized assumptions, and so it is difficult to 
see their generality because each text gives only a partial picture. It is not widely known, for 
example, that the electric and magnetic 'Hertz vectors' (which are often introduced alone, 
with their partner gauged away) transform as a skew-symmetric tensor under the Lorentz 
group, are compatible with external currents, do not require the Lorenz 13 condition, and 
have a very large gauge group containing that of the 4- vector potential. 

A comprehensive theory of Hertz potentials and their gauge freedom was developed beauti- 
fully by Nisbet [N55] lN57] and reformulated in spacetime tensor form by McCrea |M57j . (See 
also Kannenberg |Kan87j .) Nisbet's formulation, translated to the language of differential 
forms [AM T881 IT96j (which bring substantial simplification and clarity into the picture), 
was used to construct electromagnetic wavelets in [Kt)2j. We summarize the main results 
below. Please refer to that paper for details. 

Hertz potentials form a skew-symmetric tensor like the electromagnetic field itself. They are 
given in a reference frame by the electric and magnetic Hertz vectors Z e , Z m , which will be 
used here in the self-dual combination 

Z(x) = Z m (x)-iZ e (x). (132) 

These vectors are generated by electric and magnetic polarization densities P e , P m , again 
represented in the self-dual form 

P(x) = P m (x)-iP e (x). (133) 



13 



Apparently due to L.V. Lorenz and not H.A. Lorentz; see |PR84I IB 
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(P m is called the magnetization and usually denoted by M.) The two fields are connected 
by the wave equation 

DZ(x) = -P(x), D = A-dl (134) 
The electromagnetic field will be presented in the anti-self dual combination 

F(x)=D(x)+iB(x), (135) 

and Maxwell's equations become 

id t F = V x F - iV x P - i J, V ■ F = p, (136) 

where the constitutive relations (in SI units) 

D = E + P e , B = H + P m (137) 

were taken into account. We will need only polarization sources, and therefore assume 
vanishing external charge- current density: J = and p = 0. Then F is generated by Z 
according to |BW751 p. 80] 

F(x) = iCZ(x) (138) 

where £ is the operator 

CZ = V x (V x Z) +id t V x Z. (139) 

External charge-current densities can be included in (jl34j) through stream potentials |N55j . 

As in the scalar case, we first construct sourceless EM wavelets. These will split into advanced 
and retarded parts, which are then the vectorial counterparts of the scalar pulsed-beam 
wavelets. Thus begin with P = 0, so that 

UZ{x) = Z(x) = I dk e lk - x C(k). (140) 

Jc 

Since k ■ x = k ■ x — ut, ()138|) gives 

F(x) =i I dk e ikx {-k x(fex() + mok x Q ■ (141) 
Jc 

To check that this indeed solves Maxwell's equations, note that ()136|) reduces to 

id t F = VxF ^uf = ikxf or S(Jfe) f(k) = f(k), (142) 
where S(k) : C 3 — >• C 3 is the 3x3 matrix function on C defined by 

S(k)v = inxv, n(k)=k/u, n 2 = 1 (143) 



S v = v — n(n ■ v) 



§ 3 = S. 



(144) 
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Equation ()141j) therefore reads 

F(x) =2i [dke ik - x u 2 F{k)C{k), where P(Jfe) = ~ ( § 2 + S) . (145) 
Jc 2 

By ()144j) . S(fc) has the nondegenerate spectrum {1,0,-1}, and ()142|) requires / to have 
eignevalue l. 14 But 

§P = P = P* = p 2 , (146) 

so F(k) is precisely the orthogonal projection to the eigenspace with eignevalue 1. This 
shows how Hertz potentials work in Fourier space for free fields. 

By (I145j) . F has coefficient function 

/ = 2^ 2 PC=P/. (147) 
Next, extend Z(x) and F(x) to T with the analytic-signal transform, 

Z{z) =u [ dk ®{uu) e ik - z C{k) (148) 
Jc 

F{z)=u [ dke(cuu)e lk - z F(k)f(k). (149) 
Jc 

The positive and negative-frequency parts of / also have positive and negative helicities 
[K94j, so the restrictions of F(z) to T + and 71 are positive and negative- helicity solutions. 

To construct the wavelets we need a Hilbert space of solutions. The inner product is uniquely 
determined in Fourier space up to a constant by Lorentz invariance to be 

(F 1 \F 2 )= ( - 2 f 1 {k)*f 2 {k) = A [ rf^ 2 dPC 2 . 
Jc J c 

Denote the Hilbert space of all solutions with finite norm by 

H = {F : \\F\\ 2 = (F\F) < oo}. (150) 

The wavelets will be dyadics, and to streamline the notation, we rewrite the inner product 
by thinking of \Fi) as an (infinite-dimensional) column vector and (F\\ as its adjoint 
row vector with respect to the above inner product: 

= F* : H —> C, \F 2 ) = F 2 : C ^ H, (F l \F 2 ) = F\F 2 , (151) 

where F 2 : C — ► 7i is the map of scalar multiplication c i— > cF 2 . This is simply Dirac 
notation without bras and kets, called star notation in |K94j . but somewhat more flexible 
because now F{ can be replaced by an operator. 

Note that the measure dk/u 2 = dk/2\k\ 3 is invariant under scaling. In fact, an equivalent in- 
ner product has been shown to be invariant under the conformal group C of Minkowski space 

14 §/ = / means F is anti-selfdual, Sf = — f means F is selfdual, and Sf = means / is longitudinal 
and F = 0. 
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Gr64 . Therefore the Hilbert space of anti-selfdual solutions carries a unitary representation 
ofC. 



Returning to (jl49j) . define the matrix- valued function 

W z (k) = u &(uju) e~ ifc ' 2 * u 2 F(k), (152) 

represented in spacetime by 

W z (x') = u [ dk e{uu) e iHx '- z,) cu 2 P(jfc) = W(x' - z*). (153) 



c 



This is matrix- valued solution of Maxwell's equations (i.e., every column is a solution) 
designed so that its vector- valued inner product with any solution F e TC is F(z): 

r dk - 

W* F = — W z (k)*f(k) = F(z). (154) 
Jc 10 

This is a vector form of the evaluation maps (J3U|) used to define the scalar relativistic coherent 
states e z . 

The W z s are the the sourceless electromagnetic wavelets. They span a reproducing kernel 
Hilbert space with a matrix kernel 

K{z',z*)=W* z ,W z = / \ eMe(^ M )e ifc ' (2 '- 2 * ) u; 4 P(fc) 

Jc ^ 

= Q(-y'-y)W(z'-z*), z',zeT, (155) 

where the factor Q(—y' ■ y) enforces the orthogonality of wavelets parameterized by the 
forward and backward tubes and the holomorphic matrix function 

W(z)= / dk e(-k-y)e ik - z to 2 F(k) = [ dk e ik ' z u 2 P(jfe), z e T± (156) 

generates the entire wavelet family by translations: 

W z (x') = uW(x' - z*). (157) 

We now compute the 'mother wavelet' W(z) explicitly. Applying it to a vector p £ C 3 gives, 
by the same reasoning as in ()141j) and ()145j) . 

2W(z)p = 2 [ dk e ikz u 2 F(k)p= -C[G 4 (z)p], (158) 
Jc s 

where, according to (j4*Uj) and (jMj) . 

G 4 (z) = - [ dk Q(-k ■ y) e lk ' z = , 2 e T (159) 



c 
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is the extension (jSljl of the Euclidean potential G4(x E ). Since p E C 3 is arbitrary, 
determines W(z). 

Thus we have made contact with the scalar theory. We already have a good understanding 
of G±(z) and its relation to causality. By (jfi7j) . we have a splitting 



G i (z)p = iD-(z)p-iD + (z)p, D ± (z) 



1 



8«7r 2 r(r =F f) 

It is natural to define retarded and advanced Hertzian dipole pulsed-beam potentials 



(160) 



Z±{z)=D ± (z)p 



(161) 



whose polarization sources are 

P p (z) = -nZ^(z) = —pDD ± (z) =p5 3 ,i(z). 
Thus p is interpreted as a combination of magnetic and electric dipole moments 

P = Pm ~ We 



(162) 



and Zp(x—iy) are interpreted as emitted by a dipole disturbance initiated at iy and observed 



at x. By 



the associated pulsed-beam fields are 



F${z)=iCZ}{z) 



Thus ()160j) induces a causal splitting of W(z), 

W(z)p= W + (z)p- W-(z)p, 



(163) 



(164) 



where W + (z), interpreted as retarded and advanced electromagnetic pulsed-beam 
propagators, give F^ directly from p by 



2W ± (z)p = F^(z). 



(165) 



Since D ± {z) is holomorphic in T outside the world tube t> swept out by the source disk T>, 
so are Z^ , F^ and W*. The dipoles are spread over the source disk V and modulated in 
time, as specified by the distribution <5 3j i in (j^Hj) . Note that since is complex, the electric 
and magnetic dipoles become thoroughly mixed while being 'translated' from the origin to 
iy. 

The Fourier transforms in x of W ± (x — iy) can be computed from those of the scalar sources. 
For uj > 0, (ITHll and (fTT?H with g(u, u) = C(u, u) give 



C(u, u) 
2/i 



ifj,a 



-ifia 



+ 



H + Q + il ji + 0-il 



P 



(166) 
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The EM wavelets (with u > 0) are now obtained in Fourier space by taking the transform 
of (|163j) and remembering that § + S 2 7^ 2 P since we are not on-shell: 



2 W(k, y)p = + 6 \ {ik x(kxp)+ujkxp}. 



(167) 



Remarks 



• All quantities are real in spite of the the fact that we are dealing with holomorphic functions 
in complex spacetime! For example, 

D±(x - iy) = Re F£(x - iy) (168) 
Bp (x - iy) = Im F^(x - iy) 

define a real electromagnetic field in Minkowski space, for any given imaginary source point 
iy. This justifies our motto, real physics in complex spacetime. 

• The holomorphy, rooted in the cone structure of relativistic wave equations as explained in 
the introduction, serves to organize the equations by pairing dual or 'harmonically conjugate' 
fields. In the case of Maxwell fields, the duality is between magnetic and electric entities. 
This pairing survives even the introduction of sources in that local sources introduce only 
local singularities. This makes it possible to study the singular sources in terms of the 
boundary values of the fields, as we have done. 

• Note that just as holomorphy pairs electric and magnetic dependent variables, it paired 
positions and momenta as independent variables in the relativistic coherent-state represen- 
tations of massive fields. 

• It is easily shown that 

P(ife)*=P(Jfe) W(z)*=W(z*) => W z (x')* = Wf(x'), (169) 

therefore we need only consider z E T + . 

• Note that 

G{z) = s~ 2 G{z/s), s^O. (170) 
If s is independent of x (but possibly depends on y), then 

W(z) = s" 4 W(z/s). (171) 
Taking s to be the resolution parameter discussed below (J33|) . 

S=V ^=A>0, (172) 

(fT5%|) and (fTBTj) show that 

w z (x') = ±a- 4 w (^r~) = ±A ~ 4 w {^jt ~ ^) ' zeT± (173) 
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where 



y='jeV ± , f = -l. 

Thus all the wavelets are obtained from W(x — iy) with y on the hyperboloid y 2 = — 1. 
Using Lorentz invariance, y can be further restricted to y = (0, ±i), which may be further 
reduced to y = (0,i) by (jl69|) . In this way, the entire family of EM wavelets is obtained 
from a single matrix function. Furthermore, the columns and rows of W are constrained by 
()152j) since P is the projection matrix to a one-dimensional subspace. 

• On the other hand, because the splitting G4 = iD + — iD~ depends on the inertial frame 
where the branch cut is taken (the rest frame of V), we cannot apply Lorentz covariance to 
the pulsed beams W*(x'). But D ± (z) are still positive- homogeneous, 

D ± (z) = s~ 2 D ± (z/s), s>0, (174) 

therefore ()173|) still holds for W*. All these wavelets can therefore be obtained from those 
on the hyperboloid y 2 = — 1 in T + , whose elements y represent the normalized 4- velocity of 
the source disk. 



• There exist many equivalent resolutions of unity in 7i [K94J, obtained by integrating over 
various parameter sets V C T with appropriate measures d/i-p: 

[ dfi v (z)W z W* z = I n . (175) 
Jv 

This is a 'completeness relation' dual to the '(non)-orthogonality' relation ()155|) . One natural 
subset for a resolution is the Euclidean spacetime as in (@U), where all the wavelets are 
spherical and are parameterized by their center x and scale u. 

• Each resolution gives a representation of EM fields as superpositions of wavelets, 

F{x) = [ d(i P (z)W z (x?)W z F = [ du. v (z)W z (x')F(z), (176) 
Jv Jv 

with the AST F(z) restricted to V as the "wavelet transform." 

• Applying a conformal transformation to any resolution of unity gives another one. Since 
C ~ SU (2, 2) acts on T by matrix-valued Mobius transformations, the new parameter space 
is simply the transform of V. 

• The resolutions of unity considered in |K94j were all continuous, but they can be discretized, 



provided the sampling is sufficiently fine, to give frames of EM wavelets. 
• Applying the AST to ()176|) gives 



F{z')= [ dfx v (z) W* z ,W z F(z)= [ dfx v (z)K(z',z*)F(z), (177) 
Jv Jv 



'V JV 

which explains the term "reproducing kernel." 
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Combining ()164j) and ()176|) gives 

F(x') = F + (x) -F'(x') 
FV)= / d^(^W z ± (x / )F( Z ), (178) 



which is interpreted as a resolution of the sourceless field F into retarded and advanced fields 
generated by complex sources distributed over V (more precisely, on disks V> z parameterized 
by z E V). 

A great many (most!) aspects of EM wavelets, both theoretical and practical, remain un- 
explored. As previously mentioned, an exciting possibility is that the pulsed-beam wavelets 
may be realized by simulating their sources. I hope to report on the continuation of 
this research in the near future. 



10 Appendix: proofs of the source theorems 

We begin with some preliminaries concerning the oblate spheroidal coordinates associated 
with the complex distance in C 3 (for C ra , see jKOOj ). 



f(x — iy) = yjx — iy) 2 = p — iq, p > 0. 
Let V be the gradient and A the Laplacian with respect to x, for given y^O. Then 

f 2 = z 2 =>■ rVr — z, (Vf) 2 = 1. 

The unnormalized OS basis Vp, Vq is given by 

z r*z px + qy py — qx 
Vf = — = Vp = — z~z — > Vg = — — . 

ry ry^ ry ly^ ry ry^f- ry 

Its normalization and orthogonality follow from 



Vf • Vf = 1 (Vp) 2 - (Vq) 2 = 1, Vp ■ Vq = 

\z\ 2 n r r 2 4- n 2 

Vf* • Vf = ±Lr (Vp) 2 + (Vq)' 



which give 



iy t- iy iy -f- iy 



tyf-iy lyf-ry 

Taking the divergence of f Vf = z gives 

2 2p Iq 

Af = — , hence Ap = -r~z , Ag = —-^-3 . (179) 

ry ry^ ry ty^-iy 

To compute volume integrals in the oblate spheroidal coordinates, recall (|5fi|) that 

a£ = pq, a 2 p 2 = (p 2 + a 2 )(a 2 -q 2 ). (180) 
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Therefore, using differential forms (see |AMT88j . for example), 

pdp A d£ = i d(p 2 ) A d£ = d [(p 2 + a 2 ) (a 2 - q 2 )} A d(pq) (181) 

= a -1 (pdp — qdq) A (pdq + gdp) = a _1 |f| 2 dp A dg, (182) 

where dp A dg denotes the antisymmetric exterior product of differential forms. Therefore 
the volume measure in OS coordinates is 

dx = a~ l \f\ 2 dpdqdcj). (183) 



Proof of Theorem [IJ. 

Using the shorthand 



W E (z) = Q(p-ea 



g(r-r) 



6 



9 



we have, taking into account the above relations: 



(Vp) 2 = N, 



therefore 



-2nvw £ = -5 —vp + e — vf + e vf 

2r 2r 2r 2 

-2nAW £ = -5' —„N + 5 9 -N + 5%N -5-^Ap-B^. 

2r r r 2r 2r 



2ttS £ = -2nnW £ = -5'-^ z N + 5^N + 5-^N-5-^zAp. 

2r r r 2r 



Multiplying through by \r\ and letting a = p + a , 



2n\r\ S £ = — o — + + — — o — 

ry ry £i ry 



2f 



Since the only singularities are in x, no smearing needed in t and S £ acts on a test function 
f(x) by 



(S e ,f) = J dxS £ (x-ty,r)f(x) 



oo pa 



dp / dq 



_ 5 ,<?_9L + 5 °_lL + 5 <?_9f_ 5 egl 

ry <y " <y 



where f(p,q) is the mean (|91|) of f(p,q,4>) over 0. Integrating the first term by parts in p 
and simplifying gives 



S £ J] 



e 2 + a 2 ra 



2a 



dq 



g'f gf gf P 

_ "T _„ T 



;i84) 
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But g(r — r) is analytic, hence 

g '(r -r) = -g p 
Integrating the first term by parts in q gives 

~9f 



S £ J) 



a a 
~2a 



IT 



a*a f a dq ° 
a /_„ r 



(185) 



p=e 



where 



2/r fp~^~ ifq ■ 

This proves ()90p. To prove (|92|) . simply apply the expression on the right to a test function 
and integrate using (|183j) . ■ 

Note: Without the last integration by parts, we would be left with (jl84j) . which defines 
a perfectly good distribution supported on E ea without any point sources, given without 
smearing by 

S(p-ea) | a , 2 f g' , g q 



SJz) 



An 



;i86) 



However, when the limit e — > is taken in (jl84j) . the term gf/r 2 diverges as q~ 2 . The 
integral must therefore be regularized by subtracting and adding a Taylor approximation to 
the numerator gf. The resulting limit (S,f) is a nonlocal distribution, i.e., it cannot be 
written in unsmeared form like (|97j). Therefore, the last integration by parts is needed to 
obtain the local expression for S(z). 

Proof of Theorem \% 

Since the test function is continuous and f = ±iq denote the same point on T> (regarded as 
being in the upper and lower layer), we have 

/(-*?) = /(«?)■ 



Furthermore, (|56|) gives 



a 



p — 9 i 2 P 



p 2 + a 2 p 

qd p + 

a z p a 



;i87) 



2a ^ 

Of 



a — it; 



d p + id[. . 



r p 

This shows that both partials are antisymmetric in q on T>, with 

fpin) = -fp(-iq) = -ft 

a 

fq( i( l) = -fq(- i( l) = ~-fp 
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Inserting this into (j!85|) . taking limit e — > 0, and observing that 

r(z) = Ha =^ x = 0, hence f(±ia) = /(0) 



gives 



g(r + ia) g(r - ia) 



f(0) + a f ^- g ( T + iq)fr(-iq) 



ra da 

g(r,a)f(0) + 2a I -1 g( T ,q)f f (-iq) 



= g(r,a)f(0) 
Changing the integration variable to p gives 



— 9(t, q) { af p + ipfe 
o P 



(SJ}= g(r, o)/(0) + f ^ sf(r, g) W p + z>£ 

Jo Q 1 

where g = + \Ja 2 — p 2 . This proves (|9l)j). Again, (|9*7j) is proved by applying the right side to 
a test function and integrating. ■ 

Proof of Corollary [IJ. 

For g (t) = 5(t), g{r) = l/2niT and 



^ ) = 4™ 

IT 



1 1 

+ 



2nz 2 



t + iq r — iq 
z G V. 



IT 



2n(r 2 + q 2 ) 



Since W(z) = iD + (z) in this case, that proves 
For g (t) = 2, g(r) =u = g(r, q) by (JHZD and 



W(z) 
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■mG 3 (z) S(z) = uS 3 (z), 



proving (|9*9*)l . Finally, for g (X) = e g(r) = u9(wji) e K " T and 



g(<r,g) = -e(um) e — ( e ^ + e -^) 
= h0(wh) e~* wr cosh(cug). ■ 



Proof of Theorem EH. On i? £a we have 

6(g) = eg, p(g) = 



|?7| y a? — q 2 



M Po(g)- 

We will need the temporal Fourier transform of g(r — f), 

dt e^^(t -iu-r) = g{u, u) e iujf = g{u, u) e i£uja &" 1 = A 
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since the integration contour can be moved by t — > t + f without crossing the discontinuity 
of g(r — r) across the real axis. (Note that R + r — R — iq, so the deformation is bounded 
uniformly by q 2 < a 2 .) For the test function in (J9Uj) choose the plane wave 

Thus, using ffTSTI) . 

/(f) = /(p,0,£) = J (M^ = e-^Uhp) = e-^J Q (hp) 



r { p J 

Inserting this into (J9U)) and taking the temporal transform gives 
S £ (k,y) = A(u,u) Vt £ (k,y) 

n „ , \r]\ 2 a(e^ a e~^ a 1 I77I 2 T , WD f ha-iehq T/1 , , T . 



2i { rja r]*a J 2 J_ a y p 

Recalling that h £ = \rj\h and thus hp = h £ p Q , this simplifies to 

^e{k, y) = cosh(u; e a) — is sinh(cj e a) — I x + iel<i + 77*77/3 



where 



_ h £ a f a dq , . , /" a o?g .. . . . 

h = — — e Eq Ji{h £ p ) =h £ a — cos{iu £ q) Ji{h £ p ) 
2 i-a Po 7o Po 

h= he f a dq ^ = a -ig^ h 

1 J -a PO 

I f a f a 

- / dq J {h £ p ) = I \ dq cos{iu £ q)J {h £ p ). 

^ J~a JO 



h 

The first two integrals can be evaluated by letting 



g = acos7, 2tp = p £ a = a^J h 2 — uj 2 , 2% = iu £ a, ip 2 — x 2 = h 2 a 2 /4 . 
Then, using |GR651 p. 742 (6.688-1)], 

/•tt/2 

Ii = h £ a I d'y cos(iuj £ cos 7) J\{h £ a sin 7) 

= (h £ an/2)J 1/2 (ip + x)JiM ~ 
= h £ a^/i/j 2 - x 2 jo{^ + x)jo(& ~ x) 

sin(^ + x) sin(V> - x) 



v 7 ^ 2 - X 2 
2 sin 2 ^ cos 2 x — 2 cos 2 ^ sin 2 % 

2 cos 2 x — 2 cos 2 "0 
cos 2% — cos 2ip. 
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Thus 



Ii = cosh(u; e a) — cos(/i £ a) (189) 

h = = sinh(^ e a) - (w e /fi E ) sin(// e a), 

a ou £ 

where we have used 

The third integral is obtained using |GR65t p. 737 (6.677-6)]: 

J 3 = {I /lie) sin(/i e a). (190) 

Inserting these into (|188|) gives some 'miraculous' cancellations resulting in (j!05jl . The details 
are given Section [7| ■ 

11 Conclusions 

Analytic continuations to complex time and complex spacetime abound in physics, although 
the terminology of 'Wick rotations' is, in my opinion, sometimes used too casually, without 
any mathematical justification or even any basis for justification ('not even wrong'). There 
have been times while reading papers (or even books) on string theory, for example, when was 
unable to tell whether the author was working in a Euclidean or Lorentzian signature. But 
even when justified, the extensions are usually regarded as 'mathematical methods' without 
any particular physical significance. Here is a non-exhaustive list of examples known to me. 

• In the correspondence between quantum field theory and statistical mechanics, the imag- 
inary time (more precisely, its period) is related to the reciprocal temperature. But this 
is regarded as an analogy between the two theories, albeit a precise and very useful one. 
To make it more than analogy one might, for example, interpret the complex time as a 
combination of evolution and thermal parameters for a system in a local equilibrium state, 
something like the complex combination of the (also incompatible) position and momentum 
observables occurring in coherent-state representations. (These need not be eigenstates of 
a corresponding combinations of operators, as they are in the Bargmann-Segal representa- 
tion. For example, the relativistic coherent states e z (|29|) do not depend on the existence of 
'covariant spacetime operators,' which do not in fact exist within the usual framework. 

• In Wightman field theory, n-point functions are extended to tube domains in their difference 
variables and powerful methods of complex analysis are used to prove theorems like PCT 
and the connection between spin and statistics about the original fields in real spacetime 

SW64]. There is no attempt to interpret the complex coordinates z = x — iy, although 
the interpretation of y as (proportional to) an expected energy-momentum in relativistic 
coherent states, proved for free fields in |K77| IK78| IK87j , extends to general axiomatic fields 

KM Section 5.3]. 
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• In constructive quantum field theory j(XJ87 . the Euclidean region is used to correlate the 
n-point functions of a given theory by rigorous (Feynman-Kac) path integral methods. Then 
they are continued back to real spacetime and used to construct interacting fields. Again 
there is no attempt to interpret complex spacetime because the quantized field exists only 
in the Minkowskian region while the random field exists only in the Euclidean region. In 
between, there are only n-point (Wightman) functions. 

• There have been various efforts to represent spacetime as a Shilov boundary of a complex 
domain (see |G01| and references therein), but I am not aware of any claiming to do physics 
directly inside these domains. 

• Complex spacetime plays a prominent role in twistor theory |PR86t IP87j and the theories 
of Heaven or H-spaces |HJNPT78l BFP80J, but again no direct interpretation is generally 
given to the complex coordinates 15 

To the best of my knowledge, the only examples (aside from the relativistic coherent states 
and physical wavelets covered here) where complex spacetime coordinates are given a direct 
physical significance have appeared in the works of Newman et al. |N73[ INW74j . who have 
proved the following very intriguing result. 

Consider an isolated classical relativistic system in flat spacetime with positive total 'mass' 
(i.e., m 2 = Pq — P 2 > 0, where P M is the total energy-momentum). The total angular 
momentum splits into orbital and spin parts L + s, and L is made to vanish by translating 
to the center of mass. Similarly, if the system has total charge e ^ and a magnetic moment 
H, then its dipole tensor is reduced to /x by translating to the center of charge. However: 

• With a further imaginary translation by is/mc, the spin can be made to vanish. Thus 
spin may be identified with an imaginary center of mass. 

• With an imaginary translation by i/x/e, the magnetic moment can be made to vanish. 
Thus magnetic moment may be identified with an imaginary center of charge. 

• If the centers of mass and charge coincide, then the spin and magnetic moment can be 
transformed away simultaneously by an imaginary translation. The necessary and sufficient 
condition for that is that the gyromagnetic ratio of the system have the Dirac value: 

ix = (e/mc)s. 

In the massless case, the world lines with complex center of mass are replaced by a totally 
null complex plane if the spin (in real Minkowski space) is nonzero. 

This idea, although proved in flat spacetime, was inspired by the Kerr-Newman solution to 
the Einstein equation |N65j . which is the universal model for spinning, charged black holes. 
It was discovered by performing a somewhat mysterious complex coordinate transformation 
on the spherically symmetric solution with mass and charge (Reissner-Nordstrom) which is, 
roughly, a general-relativistic version of extending the Newtonian potential from IR 3 to C 3 , 
i.e., 

G*(x) = -— -> G«(z) 



Anr ' Airr 



15 



E T Newman, private communication. 
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The Kerr-Newman solution was soon realized to have the Dirac gyromagnetic ratio. Recently, 
an old debate was re-ignited with A. Trautman whether the Dirac value necessarily depended 
on the nonlinear character of the equations. Newman settled the question by showing that 
the Dirac ratio was obtained as well for the linearized solution |N02| . In the related work 

KOlaj, the charge-current distribution for a (real, static) electromagnetic field defined as in 

N73 by a holomorphic Coulomb potential 

E(x - ly) - iH {x - ly) = -V = (191) 

was computed and shown to represent a rigidly spinning disk (the branch cut T> of f) with 
radius a = \y\ and angular velocity 

u> = cy/a 2 , 

so that the rim moves at the speed of light. This is consistent with the fact that (jl91|) 
represents the electromagnetic part of the linearized Kerr-Newman black hole. 
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